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SUMMARY 


A program was conducted under NASA Contract NAS3-20292 to investigate techniques 
for simulating and evaluating the dynamic characteristics of a mixed flow turbofan engine. 

A digital computer was used to simulate the engine which was modeled using one dimen- 
sional forms of the dynamic mass, momentum, and energy equations and the dynamic 
response to planar disturbances was evaluated. The TF30 engine was simulated to allow 
comparison of dynamic characteristics to results obtained from testing of the TF30 
engine at the NASA-Lewis Research Center. During these tests, the engine was subjected 
to planar inlet pressure disturbances which were induced by the NASA-Lewis developed 
air Jet system. 

The simulation techniques investigated during this program concentrated on improving 
the dynamic characteristics of the engine simulation through modification of the compres 
sion system model. Improvement of the simulation dynamic characteristics was approached 
in three steps; First, a baseline engine simulation using overall component performance 
characteristics and dynamic elements sized to the lengths of the major compression 
system components was assembled, and the simulation dynamic characteristics were 
established. During the process of establishing the simulation dynamic characteristics, 
the importance of the dynamic energy equation was assessed. Second, a separate simulation 
of the compression system was assembled using blade and vane row performance character- 
istics. Six dynamic configurations of the compression system were formulated by combining 
the row performance characteristics to form dynamic elements composed of compressor 
stages or stage groups. The dynamic characteristics and computer run time were established 
and an optimum dynamic configuration was selected on the basis of comparison to the 
dynamic characteristics and computer run time associated with a stage-by-stage simulation. 
Computer run time improvements associated with algorithms for prediction of the program 
iterative solutions were also investigated using the optimum dynamic configuration of 
the compression system. Third, the engine simulation was modified by incorporating the 
optimum compression system. The dynamic characteristics of the modified engine simula- 
tion were then evaluated for improvements associated with modifications to the compres- 
sion system model. 

Three tests were employed to evaluate the dynamic characteristics of the simulations. These 
were tests of the simulation frequency response (gain and phase versus frequency), transient 
response (pressure versus time), and limits of system stability. Each of these were used to 
evaluate cue simulation versus another and/or to evaluate the simulations versus test results. 

Three methods were used to calculate simulation frequency response and consisted of 1) 
transfer function, 2) step response, and 3) discrete frequency methods. The transfer function 
method was different from the other two in the respect that state variables were used to de- 
scribe the simulation by a high order system of linear differential equations. Solution of the 
differential equations established the transfer functions describing the frequency response 
characteristics. The step response and discrete frequency methods employed analysis of the 
simulation input and output transients to arrive at frequency response characteristics. 



The significant results obtained in this study included the following: 

• The reduced order analysis of thci compression system simulation showed that at 
100 the dynamic characteristic of a stage*by-stage simulation could be 
obtained with a more simplified dynamic configuration while realizing substan* 
tial savings in computer execution time. The analysis showed that the dynamic 
elements should be uniformly sized to achieve the best frequency response. 

• The engine simulation with the modification to the compression system showed 
improved frequency response characteristics above 50 Hz. 

• The transfer function provided better resolution of the high frequency response 
of the simulation than the transient methods. 

• Coupling of the fan bypass and core dynamics has a potentially significant effect 
on core stream transient operation. 

• Stability limits of the modified engine simulation were twice that shown by en- 
gine test results. 

The resulto of this program highlight the need to further develop the engine simulation to 
achieve closer agreement with actual engine dynamic characteristics. 
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INTRODUCTION 


Gas turbine manufacturers must produce engines which operate stably during normal use. 

During normal operation, the engine is exposed to time variant operating conditions such as 
gas generator power transients, augmentor transients, and inlet distortion which can poten- 
tially produce compression system surge/stall. The equivalent dynamic frequency of these 
time variant operating conditions can range from less than one to thousands of Hertz. Engine 
response to the time variant conditions is dependent on the coupled dynamic characteristics 
of the assembled system of engine components which in turn affects stability. 

To maximize the stability of a developed engine, it is necessary to assess engine stability 
characteristics during the design stage. To accomplish this, reliable analytical capability is 
required. Computer programs which simulate the low frequency dynamics of engine power 
transients are being used successfully in support of the design process. In contrast, dynamic 
simulation of engine operation in a high frequency environment such as time variant inlet 
distortion, is not as reliable. As a result, the engine design criteria for stability are primarily 
derived from past experience and established within the framework of uncoupled component 
characteristics. The earliest that stability of the engine as a dynamic system can be evaluat- 
ed is during the development stage. Since the cost of the instrumentation and tests is often 
high, the testing is usually of limited scope, and as a result, the impact on the engine design 
is minimal. 

Simulation of engine dynamics can be achieved using analog, digital or hybrid (analog plus 
digital) computers. The analog approach provides continuous transient results with a mini- 
mum of implementation effort, while requiring large, complex and expensive computing sys- 
tems whose reliability and accuracy decreases with the size of the program. Digital computa- 
tion provides the capability to model large systems with a high degree of precision at the 
expense of computer run time and increased implementation complexity. Hybrid systems 
have been formulated to simultaneously exploit the advantages of analog and digital compu- 
ters. In recent years, improved simulation techniques coupled with faster digital processors 
have made digital simulation of high frequency engine transients the most practical alterna- 
tive available. 

During the yean 1970 - 1974, under contract to the Air Force, Pratt & Whitney Aircraft 
assembled a digital dynamic simulation of the TF30 engine and investigated its dynamic 
characteristics. A» reported in Vol. II of AFAPL-TR-74-107, the simulation used only the 
mass equation in dynamic form and modeled the compression system using one dynamic 
element for each compression system component. Evaluation of the simulation dynamics 
revealed that modifications to the compression system simulation were required to improve 
the core stream dynamic characteristics. Stage-by-stage modeling of the compression sys- 
tem was suggested as approach for improving the core dynamic characteristics, but was 
not tried in the engine simulation because of limitations imposed at that time by computer 
size and run time. 

The work discussed in this report has used the improved computational capabilities of the 
digital computer to further investigate dynamic engine simulations. The objective of the work 
was to assemble a digital computer simulation of a mixed flow, twin spool turbofan engine and 
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evaluate the dynamic diaracteristics of the engine. The goal of the woric was to improve the 
dynamic characteristic of the simulation to disturbance ftequendes of at least 100 Hz by 
modifying the compression system simulation. The engine was modeled using one dimensional 
forms of the dynamic mass momentum and energy equations and response to planar distur- 
bances was evaluated. The TF30 engine was simulated so that dynamic characteristics could 
be evaluated against results obtained from testing of the TF30 engine at the NASA Lewis 
Research Center. During these tests, the engine was subjected to planar inlet pressure distur- 
bances induced by the air jet system developed by the Lewis Research Center. 

A three step approach was used in working toward the goal of improving dynamic charac- 
teristics of the engine simulation. First, a baseline engine simulation was assembled using 
dynamic elements sized on the order of the compression system compcment lengths and the 
dynamic characteristics established. Second, a separate ccnnpression system simulation was 
assembled by adding dynamics to the compression system simulation used for the distor- 
tion modeling of the TF30-P-3 compression system reported in NASA CR-135 1 24. The dy- 
namic element sizing of the simulation was varied by combining compressor row performance 
characteristics to establish a configuration which produced dynamic characteristics similar 
to a stage by stage simulation, while minimizing the computer computation time. Third, 
the engine simulation was modified by introducing the compression system simulation es- 
tablished in tlie second step. Hie dynamic characteristics of the modified engine simulation 
were then reevaluated for improvements. The analysis of simulation dynamics was performed 
in each of the three steps at engine operating conditions of approximately 8600 and 7600 
rpm corrected low spool speed (Nl/yjjJ^. 


engine simulation 


NASA LEWIS TESTS OF ENGINE DYNAMICS 


* ^ K NASA Lewis Research Center to investigate engine dy- 
A TF30-P-3 engine was tested by the NAS dynamic 

namics. The results obtamed from the ^^formed in an altitude chamber usmg 

characteristics of the engine simulation. T s ^escribed in ref. 1, which induced dyna* 

the Uwis Research Center develo^d i„tet conditions for the test produced 

mic pressure conditions at the engme ^07 to 5.45 N/CM^ abs (7.35 to 

a Reynolds number index of 0.5 at an ® P ^ N/CM^ abs (2.4 psia). As illustrated 
7.90 W The engine exhaust pressure ,„p„e inlet 

in Figure 1, high response total pressure in discharge (station 2.3). 

(Station 2.0). fan tirst s««or ^ ^L”re-t discharge (sUtion 3.0). twelfth 

A portion of the NASA test program ST 

sponse to spatially uniform (planar) pen^ic reference 1 gain and phase char- 

let. By using a frequency sweep test techmque at the 

acteristics otId‘’d!a“Gain and phase characteristics acquit^ at 

e.S:c p:.rc:— 0 ^ appri^^^^^^ 8600 and and 

ZJko> were us^ for and OP2 

modified engine simulations. These P 8 ^ investigated under the NASA 

respectively, in this report. ^ inlet pressure oscillations were induced 

test program. During these tests, dis^te or the ampUtude 

and ampUtude was increased until either engine s J ^ p ampUtude 

limits of the air jet system were reached ^es * ^ modified^ngine simula- 

at staU/surge were used for comp ^sonto the stabiUty limits 

tion at approximately 8600 N1/^0'T2.O’ 

SIMULATION SYSTEM 

Pratt A Whitney Aircraft has was to provide the cap- 

(State of the Art Performance Program). ^ (-.^-nons with simple commands to the 
ability to assemble simulations of new enpne cor^nguratiom^simp 

System. To implement the system, basic ra"omP«t«t 

major engine component i. e. comprewr. , ' characteristics of 

library. A second Ubrary was also ™ 't*' Te ^LuCnTtem was formulated 

each component of e«ry engme ™ • - ^^yityto simulate operation of any configu- 

Si2.CC— 

tions discussed in this report. 
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BaseHne Dynamic Engine Simulation 

The first step to improving engine simulation dynamic characteristics included assembly of 
a baseline engine simulation and defimtion of the dynamic characteristics. The purpose of 
this step was to establish a baseline of dynamic characteristics against which the influence 
of compression system model modifications could be assessed. The discussion which follows 
provides a comparative description of the TF30 engine and the configuration of the baseline 
engine simulation. This is followed by an evaluation of the importance of the energy equa* 
tion to simulation dynamic characteristics. This discussion is concluded with a description 
of the dynamic characteristics derived from the simulation. 

Configuration 


The TF30-P-3 is a twin spool augmented, turbofan engine which mixes bypass and core gas 
flows in the afterburner duct. The low pressure spool is composed of a three stage fan and 
rix stage low pressure compressor driven by a three stage turbine. The high pressure spool 
is composed of a seven stage high pressure compressor driven by a one stage air-cooled tur- 
bine. The combustor is a can type surrounded by a secondary chamber which uses high com- 
pressor discharge air to cool the burner cans, and dilute the fuel-air mixture. Surge bleeds are 
incorporated at the seventh and twelfth stages of the compression system. 


The dynarmc configuration assembled for the baseline engine simulation is compared to a 
side view of the engine in Figure 2. The configuration is composed of twenty dynamic ele- 
ments sized on the order of the length dimensions of the mgjor compression system compon- 
ents (i. e. fan, low compressor, high compressor). To achieve the length similarity, the high 
and low pressure turbines were combined into a single dynamic element, while the fan dis- 
charge and afterburner ducts were subdivided into a total of 1 2 elements. The subdivision 
of the afterburner duct included separation of the bypass and core gas streams. The main 
burner was represented by two parallel elements which separated the primary burning and 
secondary cooling airflows. 

The afterburner duct bypass and core flow streams were separated to accomodate the dynam- 
ic mixing calculations investigated under the Air Force Contract reported in Vol. II of AFAPL- 
TR-47-107. The calculation assumes that during transients, the inertial and viscous character- 
istics of the gas flow induce a lag in the equalization of static pressures in the bypass and core 
streams. The calculation allows a pressure imbalance to occur during transients and equalizes 
the pressures during steady state operation. During transients the pressure imbalance changes 
the flow areas in a direction to equalize the static pressures. 

Appendix A describes the methods used to calculate dynamic operation of the modified en- 
gine ^ulation assembled for the third step of this analysis. With the exception of the com- 
pression system, the comprment performance characteristics described in Appendix are the 
same for the baseline simulation; therefore only the compression system characteristics are 
described at this point of the discussion. 
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The perfoimance characteristics of the low and high compressor are defined by the curves 
illustrated in Figures 3 and 4 which show the compressor inlet corrected airflow and tempera- 
ture rise correlated with component inlet corrected speed and m. As illustrated in Figure 5, 
the m values identify lines of constant slope on a standard pressure ratio versus corrected air- 
flow map. The values of m are calculated from the equation 


PR - CpR 
m = 


( 1 ) 


WAC - 


where 

PR - pressure ratio 


WAC = corrected inlet flow 

CpR and ~ the pressure ratio and airflow 

at the intersection of all in lines 

The m line method of defining maps is used to accomodate the regions of the compressor 
map where constant speed characteristics are either vertical or horizontal. 


Tire performance of the fan bypass stream was defined in the same m line format as the 
low and high compressors. Total inlet corrected airflow, and bypass stream temperature rise 
(ATj/Tj bypass^ define curves in the format of Figures 3 and 4 respectively. 

The fan core stream pressure ratio (PRcoRE^ temperature ratio (TR^oRE^ 
scribed as illustrated in Figures 6 and 7 respectively. 


Dynamic Energy Equation 

Prior to establishing the dynamic characteristics of the baseline engine simulation, the influ- 
ence of the dynamic energy equation on simulation transient characteristics was investigated. 
The analysis was made to determine if the energy equation was required, in addition to the 
dynamic mass and momentum equations, to define the dynamic characteristics of the base- 
line simulation. To assess the role of the energy equation, two transient test cases were com- 
puted for a step increase in engine inlet pressure. In one case only the dynamic mass and mo- 
mentum equations were used; the second case used dynamic mass, momentum, and energy 
equations. 

The results of the two test cases are illustrated on Figure 8. The figure shows the inlet pres- 
sure step (Pj 2 0^ pressure and temperature at the fan duct discharge station (7.0F) 

plotted versus time. Comparison of the two cases shows that without the energy equation, 
the Pj 2 0 induced an initial decrease in the duct pressure and temperature which was 
immediately propagated over the period of the first computed time step (DT = .001 seconds). 
This result was dictated by the fan map characteristics. During the first time step of the tran- 
sient (.001 seconds), the fan pressure ratio and consequently temperature ratio were reduced. 
The reduced temperature ratio in turn lowered the fan discharge temperature because the in- 
let temperature to the engine was held constant during the transient. In the absence of the 
energy equation, the reduced fan discharge temperature was instantaneously propagated to 
station 7.0F in the fan duct. Solution of the gas state equation at station 7. OF required a re- 
duction in the pressure level in the presence of the reduced temperature. During subsequent 
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time steps, increased airflow and pressure levels induced by the step in inlet pressure were 
propagated through the fan ducts and caused the pressure and temperature rise at station 7.0F 


The case which included the energy equation showed that a 0.003 to 0.004 second time per- 
iod elapsed before the fan duct discharge pressure and temperature responded; this compared 
to an estimated propagation time of approximately 0.005 seconds based on sonic velocity. 
airOow velocity and length of the duct. The initial change in pressure and tempera^re were 
in an increasing direction as expected. This result was attributoble 

equation which transported the fan discharge temperature changes through the <*“Ct at e^ 
tSlly the airflow velocity. On the basis of this comparison, the energy equation was mcluded 
in the simulation to achieve the most realistic dynamic characteristics. 


Dynamic Characteristics 


The dynamic characteristics of the baseline simulation were established in two formats. First, 
gain and phase plots were deflned for re^onse to Pt 2.0 ^ pressures at station 2.3, 

2.3F, 3.0 and 4.0. The gain and phase characteristics were compared to similar results from 
the NASA TF30/air jet testing. Second, pressure time traces were deflned for stations 2.3, 
2.3F, 3.0 and 4.0 for a step input of Pj^O compari- 

son to the modified engine simulation. 


The frequency response, i. e. gain and phase versus frequency of the simulation was evaluated 
by three methods: fransfer functions, direct analysis of the Pt 2.0 transient response, 
and amplitude and phase analysis of response to discrete frequency Pj2.0 waves. The 
capability for all of these analytical methods was consolidated into the engine simulation and 
could be activated via input. 

The transfer function was the first method used for computing response characteristics. This 
method was implemented by combining a state variable definition of the engine with the 
routines of Seidel reported in reference 2. As described in reference 3., the state viable 
method provides a linear description of a system using first order differential equations of the 
form. 


jl(t) = Ax(t) + Bu(t) where 
X = state vector 
u = input vector 
A = system matrix 
B » input matrix 

The system output is described by 

y(t) = Cx(t) + Du(t) where (3) 

y ■ output vector 
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C = output matrix 

D = output/input direct coupling matrix 

The matrices in the above expressions consist of partial derivatives which are determined from 
the simulation by finite differences. Each state and input variable was perturbed individually 
vrith all othen held constant. 

The state variables used for the baseline simulation were the mass, momentum, and energy of 
each dynamic element, core exhaust area for dynamic mixing and fuel flow for the engine 
fuel control. The input variable was Pj 2 0 output variables were Pj 2 3 , Pj 2 3p’ 

**T3. Q and **T4.0- Formulation of these variables resulted in matrices of the following form 
for a 62nd order system. 


ail 

3xj 

B = 

’i, 

3x1 

1 

1 

’*62 

1 

i 

1 

3 u, 
1 1 

1 

1 

1 

1 

1 


1 

\ 

’*62 

3xi 

^H2 


3 U] 

3yi 

3yj 

D = 

3y| 

®,*i 

1 

1 

1 

’*62 

1 

1 

’ll, 

1 

1 

1 

>. 

ro 

.1 


1 

1 

9y4 

3x1 

’*62 
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Using the computer routines of Seidel which were furnished by the NASA Lewis Research 
Center and reported in reference 2, the matrix algebra was performed for simultaneous solu- 
tion of equations (2) and (3) to define the transfer functions-^^. The resulting gain and 
phase characteristics of the simulation were obtained from the transfer functions. 

The second method for defining simulation frequency response used the pressure versus time 
characteristics computed for the Pj 2 q step transient. Using the methods of Levy reported 
in reference 4 the gain and phase were computed directly from the time response of Pj 2 31 

**T2.3F* **T3.0 **T4.0- 
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The third method for deriving gain and phase characteristics used direct analysis of a transient 
induced with discrete frequency sinusoidal waves. As illustrated in Figure 9, direct 
measurement of the input (Px2.0^ output 0^2.3’ ^T2.3F» ^T3.0 ^T4.0^ 

plitudes and output lag were made from the pressure time traces of the transient. A half wave 
amplitude of approximately 0.25% P>p 2 q was used for input and produced smoothly defmed 
output sine waves. The simulation transient was run for a period of time which insured that 
stabilized output wave amplitudes and lags were reached. The ratio of output over input am- 
plitudes was calculated to define gain, and the lag time was converted to phase angle. 

Gain and phase response characteristics were computed for the baseline engine simulation 
at conditions OPl and OP2. The average absolute inlet and exhaust pressures were set at 
5.34 N/cm^ (7.75 psia) and 1.65 N/cmr (2.4 psia), respectively, at condition OPl to simulate 
the conditions of the NASA engine tests. At condition OP2 the inlet and exhaust pressures 
were 5.44 N/cm^ (7.89 psia) and 1.65 N/cm^ (2.4 psia). 

The response characteristics derived at OPl using the three methods discussed above are shown 
on Figures 10 through 17 and are compared to the engine test results. The gain levels were 
normalized to steady state compression ratio and plotted in units of decibels versus input 
frequency in units of Hertz. Phase was plotted in units of degrees versus frequency. Figure 
10 shows the gain response of Pj 4 q to Pj 2 q- levels derived by the transfer function 

and step response methods are shown over the range of 0. 1 to 200 Hz. Gain for the discrete 
frequency cases are shown as points at 20, 35, 54, 73, and 100 Hz. The gain levels computed 
by the discrete frequency wave and step response analyses of the induced transients show 
good agreement. Gain levels computed by the transfer function method were on the order 
of 1 to 2 db higher than the levels derived from the transients at frequencies up to 50 Hz; 
at higher frequencies the levels diverged. Comparison of the simulation gain levels to the 
engine data shows good agreement to 10 Hz, with the simulation levels signiflcantly higher 
above 10 Hz. 

The phase response of Pj^ q to Pj 2 o ** shown for the simulation and engine test data on 
Figure 1 1. The phase angles computed from induced Pj2.0 ^*‘**'sients (step and discrete fre- 
quency) formed a band of results which was up to 30® wide at frequencies up to 100 Hz. The 
transfer function phase levels were to the lower side of the band in close agreement with dis- 
crete frequency results. Comparison to test data shows good agreement (~ 1 5® A phase) up 
to 10 Hz. Above 10 Hz, the simulation and test data begin to show significant divergence. 

Figures 12 through 15 show the gain and phase response for the remaining core flow stations 
2.3 and 3.0 at condition OPl. These curves shows results similar to station 4.0. Simulation 
gain levels determined by the three analytical methods are in close agreement at frequencies 
up to 50 Hz and in close agreement with test data up to 10-20 Hz. Simulation phase charac- 
teristics show banded results to 100 Hz and significant deviation from test data at frequencies 
above 20 Hz. 

The gain characteristics of the simulation and test data at the fan discharge station 2.3F are 
illustrated in Figure 16 for condition OPl. The simulation gain characteristics computed by 
all three methods are in good agreement up to a frequency of 65 Hz. All methods predict a 
response peak at approximately 35 Hz which is close to the 36.6 Hz equivalent frequency 
estimated for downstream plus return propagation of a pressure disturbance in the fan ducts. 
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Above 65 Hz, the transfer function method computes response peaks at approximately 74, 

1 18 and 188 Hz which are close to the 36.6 Hz first, second, and fourth harmonics of 73.2, 
109.8 and 183 Hz. In contrast, the gains computed from the step and discrete frequency 
transients show only slight resolution of a peak at 74 Hz and complete attenuation of higher 
order peaks. The engine data shown in Figure 16 shows response peaks at approximately 32 
Hz and 82 Hz. Data was defined to only 1 00 Hz and it is therefore not known if other peaks 
would be resolved above this frequency. The close correspondence of the test peaks at 32 
and 82 Hz to the estimated equivalent propagation frequency of 36.6 and the first harmonic 
frequency of 73.2 Hz led to the conclusion that the data was showing a dominating influence 
of the fan duct dynamics. Comparison of the simulation results to the test data indicated 
that the transfer function method gave a better description of the actual predominant fan 
discharge dynamics than was shown by transient computations. 

The station 2.3F phase characteristics of the simulation and test data at condition OPl are 
shown in Figure 1 7. As in the case of the core stations, the simulation phase characteristics 
computed from inlet pressure transients gave a band of results which showed a spread of up 
to 35° at 100 Hz. The phase derived from transfer functions again fell along the bottom of 
the band in good agreement with the discrete frequency transients. The transfer function and 
step transient methods predicted phase peaks at approximately 25 and 62 Hz, however, peak 
attenuation was again in evidence by 62 Hz. Phase peaks were not identified from the dis- 
crete frequency results. It is expected that the phase peaks would not be resolved by addi- 
tional discrete frequency cases, but would show the smoothing demonstrated by the step 
transient results above 35 Hz. This expectation is based on the good agreement of the gain 
levels predicted by the two methods as shown on Figure 16, plus past experience which 
has shown similar results from both methods. The engine test data and simulation predictions 
showed good agreement at the 25 Hz peak. Above 40 Hz divergence of engine and simulation 
characteristics occurred with the second phase peak of the test data occurring at approxim- 
ately 90 Hz, compared to 62 Hz for the simulation. At condition OP2 simulation gain and 
phase characteristics were generated using only the transfer function and step response me- 
thods. The simulation characteristics are compared to the engine test data on Figures 18 
through 25. The results obtained at OP2 generally agreed with the findings at OPl . 

In addition to the gain and phase characteristics derived at OPl and OP2, pressure and rotor 
speed versus time plots were generated for inlet pressure and fuel flow step transients. The 
results of these cases are not discussed here because they were used to support the analysis 
of the compression system simulation and modified engine simulation. 

On the basis of the results obtained from the baseline engine simulation, the following con- 
clusions were reached at this point of analysis. 

1 . Dynamic mass, momentun and energy equations should be used in higii frequency 
engine simulations. 

2. Gain characteristics obtained by the transfer function method were representative 
of simulation transient dynamics at frequencies up to 50 Hz. 

3. Phase characteristics obtained from transfer functions were representative of sim- 
ulation transient results to 100 Hz. 
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4. The transient method is better than the discrete frequency method of deriving gain 
and phase characteristics from simulation transients because full sweep informa- 
tion is provided. 

5. Simulation core stream response characteristics were representative of the engine 
at frequencies between 0 and 10-20 Hz. 

6. Simulation fan stream response characteristics were representative of the engine 
at frequencies between 0 and 40-65 Hz. 

7. Frequency response results confutned the need to construct a better simulation 
of the core airstream. 

Modifled Compression System Simulation 

The second step toward improving the dynamic characteristics of the engine simulation was 
concentrated on modifying the compression system model. The modification was based on 
the past experience gained from the Air Force sponsored program. These investigations 
suggested that stage-by-stage modeling of the compression system should be implemented 
to improve engine simulation dynamic response in the engine core airstream. The potential 
for exceeding available computer storage and introducing excessively long run time require- 
ments was recognized as a drawback to implementing a stage-by-stage model. 

On the basis of the Air Force studies, modifications to the cennpression system model were 
performed within the framework of a separate compression system simulation. Due to the 
absence of other engine components, the system dynamics were different from the engine, 
and the data acquired from NASA engine tests could not easily be used as a basis for assess- 
ing modifications; instead, a stage-by-stage simulation was assembled and used to establish 
a basis for assessing modifications. 

Modification of the compression system simulation was approached in three steps. First, a 
stage-by-stage simulation of the fan, low and high compressors was assembled by combining 
row performance characteristics to create a dynamic element for each stage. The dynamic 
characteristics and computer run time were then defined at conditions OPl and OP2. Second , 
five reduced order simulations were assembled by combining stages to reduce the number of 
dynamic elements. The dynamic characteristics and computer run times of the reduced order 
models were compared to the stage-by-stage results to select a configuration for further eval- 
uation in the engine simulation. Third, a further reduction in computer run time was investi- 
gated by evaluating four algorithms for prediction of program iterative solutions. 

The discussion which follows presents a description of the stage-by-stage simulation and its 
dynamic characteristics. The reduced order compression system modeb are then described 
and followed by a comparative discussion of the dynamic characteristics and run time re- 
quirements of the reduced order and stage-by-stage models. The discussion is concluded with 
a presentation of the predictor algorithm investigation. 
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Stage-by Stage Simulation 


The stage-by-stage simulation was assembled by adding finite element dynamics to the 
compression system simulation used for the circumferential distortion modeling of the 
TF30-P-3 as reported in NASA CR-135124. The dynamic configuration is illustrated in 
Figure 26. Key features of the model included use of rotor and stator row performance 
characteristics, accounting of bypass ratio effects on radial streamline bending in the fan and 
airflow separation off of the flow splitter at the low compressor entrance, and use of unsteady 
aerodynamic effects on airfoil loss characteristics. Due to the absence of other engine com 
ponents downstream of the compression system, choked nozzles were used to terminate 
both the core and bypass streams of the simulation. A brief description of the implementa- 
tion of these effects is provided here and supported with further details in Appendix A. 


The steady state performance was described by characteristics of 4>, and X. The perform- 
ance characteristics were derived from rig tests of the compression system and smoothed for 
data errors. The rotor and stator row static pressure rise map was described by ^ and ^ where 


(N/v/0^ DESIGN 
5>pA (N^/5y) 


^ = 



(N/v/^T> design 

(N/v/5^) 


(5) 


( 6 ) 


X = 


ATy 

Tt 


(N/v/^j) design 
(N^/5j) 


(7) 


The total pressure required for calculation of and the static pressure in ^ were related by 





y- 1.0 

+ 

2.0 



( 8 ) 


The effects of bypass ratio variations on the flow properties of the fan and low compressor 
were accounted for in the simulation and are illustrated in Figure 27. The location of the 
streamline dividing the bypass and core airflows is shown for operating conditions of high 
and low bypass ratio. Under conditions of high bypass, the streamline is located toward the 
fan hub at the plane of the inlet guide vanes. As flow precedes through the fan toward the 
splitter, the streamline bend" radially in the direction of the fan tip and stagnates on the 
splitter. At conditions of low bypass ratio, the dividing streamline starte toward the fan tip 
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it the IGV plane and bends radially toward the hub as the spUtter U approached. When the 
bypass ratio is sufficiently low, the core flow separates from the drooped leading edge of the 
sputter and reduces the effective flow area through the core stream third stator. 

The bypass ratio effects are accounted for in the simulation in three ways. First. ^ spUtting 

streamline location in the IGV through rotor 3 are computed 

unit area thiougli the fan. Second, the rotor 3 map * values are modified 

ial streamline bending. Third, core stator 3 flow area, inlet flow angle, map * value, and sah 

pressure are modified to account for effects of flow separation at the sphtter. 

Unsteady aerodynamic effects were incorporated to account for the deviations f™™ 

steady Uft characteristics that an airfoil exhibits in a transient 

separation from an airfoU (stafi) occur at iarge incidence angles and requires a 

of time to propogate to the trailing edge. The effect on airfoil ^rfo|mance » 

lire 28. The difference in airfoil Uft coefficient translates in ^ 

loss characteristics for compressor operation as lUiistrated in Figure 29. The net effwt ™ 

compressor operation is to effectively modify the compressor map characttraOcs m a 1^ 

frequency operating environment. The unsteady effects were incorporated into the compres- 

sion system simulation as descnbed in Appendix A. 

Stage-by-Stage Dynamic Characteristics 

The dynamic characteristics of the stage by stage compression system were 

the high compressor discharge sUtion (4.0). Gain and pha« P'®“ 

for a step in inlet pressure were defined for Pj 4 q- Both the transfer function and step re- 

sponse methods were used to establish the gain and phase characteristics. 

Figures 30 and 31 show gain and phase plots for the simulation at condWon OFF ^ 
characteristics defined by the transfer function and step response metho^ “1^?^ “ 
lire 30 to agree closely up to 1 20 Hz, beyond which the traults diverge. Figure M 
good agree^t in phase to 200 Hz with no more than 20 spread !«7“" 

Figures 32 and 33 show the OP2 gain and phase charactenstics-denved by the Wo ">«'*'<« 
to be in close agreement to 200 Hz. On the basis of the results, the gam and phase chwttr- 
istics derived by the transfer function method was concluded to be represen^ve ofsOTuU- 

“rS^tsL frequencies of at least 120 Hz. The transfer function meth^ was seated 

for gain and phase evaluation of the reduced order models becaure the me ttod yp ear^ to 
give^tter resolution of fan high frequency response ch^tenstira m thebaseh ne ^ 
tiroulation. The pressure-time histories of response to the step m inlet presture are shown 
and discussed with results from the reduced order simulations. 

Description Of Reduced Order Compression System Simulations 

Five reduced order compression system simuUtions were assembled and cval«at^ ^ 
^tential of reducing ctmputcr runtime while retaining the 

stoge^iy^tage simulation. Each of the reduced order simulations ^ the row I**^"** 
chwicteristia. and bypass ratio effects Implemented in the stagc-by-stafc 
duced Older configurations were esUblishcd by combining 

characteristics into dynamic elements to form stage groups within a single finite element. 
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The configurations of the five reduced order simulations are summarized in Table I. The 
five configurations were identified as: 1 ) fan/low/high lumped, 2) ten similarly sized ele- 
ments, 3) bleed, 4) fan lumped/compressors staged, and 5) fan staged/compressors lumped. 

The configurations fell into two categories; the first 3 configurations resulted in dynamic 
elements of essentially similar length while the last two resulted in dissimilar lengths. Each 
of the dynamic configurations is illustrated in Figures 34 through 38, where the dynamic 
stage groupings are shown as a function of the cumulative length of the compression system. 
Dynamic element sizing was necessarily constrained by the dimensions of each compressor 
stage. 

Reduced Order Dynamic Characteristics 

The dynamic characteristics of each reduced order simulation were evaluated on the basis of 
gain and phase plots derived by the transfer function method and pressure-time response to 
a step in inlet pressure (Pj 2 o^- Dynamics of the high pressure compressor discharge pressure 
(Pt 4 q) were evaluated on each bnsis. The computer execution time for each case was also 
established to assess the trade-off between dynamic characteristics and computer run time. 

The gain and phase characteristics derived from the five reduced order simulations are com- 
pared to the stage-by-stage results on Figures 39 through 42. Results at condition OPl are 
shown on Figures 39 and 40; results at OP2 are shown on Figures 41 and 42. In general, 
these figures show that as the order of the simulation is reduced, the gain and phase angle at 
constant frequency is decreased. 

Since the frequency responses of all the reduced order models are similar, it is difficult to 
determine which one should be used. In order to obtain a qualitative measure of the dynam- 
ics of the reduced order models, the gain and phase curves of Figures 39 through 42 were 
cross-plotted at 100 Hz. Figure 43 shows the deviation from the stage-by-stage simulation 
gain level in units of decibels plotted as a function of the number of simulation dynamic ele- 
ments. Two characteristics of the plot are outstanding: first, the stage-by-stage gain level 
was best maintained when the relative sizing of dynamic elements was similar (4, 6, 10 and 
21 element configurations versus the compressor lumped and fan lumped); second, the gain 
level of the stage-by-stage simulation (21 elements) was essentially retained when the system 
was reduced to ten dynamic elements of similar relative sizing. A third observation from the 
plot notes that the c».mpressor lumped and fan lumped points for condition OPl produced 
spurious results. These two points failed to reach solutions within the prescribed program 
iteration tolerances, resulting in poorer calculation precision than the other cases. This ten- 
dency toward iteration failures appeared to correlate with systems composed of dissimilar 
dynamic element sizing. 

The crossplot of phase results at 100 Hz is shown on Figure 44. Here, the deviation from the 
stage-by-sUge simulation phase angle is plotted in degrees versus the number of dynamic ele- 
ments. Again, results were grouped according to similarity of dynamic element sizing, with 
least phase deviation (- 20®) shown by the configuration with 10 similarly sized elements, 
rhe two points which failed to reach iterative solutions with in program tolerances again 
^owed spurious results. 
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TABLE I 


Configuration 


1. Fan/low/high lumped 

2 . Ten siinilarly sized 
elements 

3. Bleed 

4 . Fan lumped/comptessor 
staged 

5. Fan staged/compressor 
lumped 


Summary of Reduced Order Compression System Simulations 
Relative 

Number of Element State Variable 


Dynamic Elements 

Length 

System Order 

Comments 

4 

Similar 

12 

Same dynamic configuration as baseline 
engine simulation 

10 

Similar 

30 

Approximately half the dynamic elements 
of the stage-by-stage simulation 

6 

Similar 

18 

Dynamic element boundaries placed at the 
component interface and bleed stations 

.15 

Dissimilar 

45 

Fan bypass and core stages each combined 
into one dynamic element 

10 

Dissimilar 

30 

Low and high compressor stages each 
combined into one dynamic element 


Simulation pressure-time response to a step in inlet pressure is illustrated on Figures 45 and 
46. Plots of high compressor discharge pressure versus time are shown for the Five reduced 
order and stage-by-stage simulations. Results at condition OPl are shown on Figure 45: 
Figure 46 shows OP2 results. On each plot, 63.2% of the pressure change at station 4.0 is 
identified to give a measure of the effective first order time constant of the simulations. 

Here also the responses of all the reduced order models are similar. In order to make a com- 
parison of the step responses results, a crossplot of the time to 63.2^^ of the compressor 
discharge pressure rise was made as shown on Figure 47. The crossplot shows once again, 
that the configurations with dissimilar element sizing gave results outside the family obtained 
with similar element sizing. Both the bleed and ten element configurations showed a min- 
imal effect on response time, with the ten element configuration showing a response time 
closest to the stage-by-stage. 

Computer Run Time 

The impact of the simulation order on computer run time is shown on Figure 48 for the in- 
let pressure step transient. The case execution time (excluding load time) is plotted versus 
the number of dynamic elements. The curve shows that by reducing the order of the model 
from 21 elements to a range of four to ten elements, significant reductions in computed 
run time can be achieved. As a minimum, the run time reduction is approximately 55% 
with the ten element configuration. Approximately 10% additional run time reduction can 
be achieved by further reducing the number of elements from ten to four. 

Figure 49 shows the run times required for computation of frequency response by the trans- 
fer function method. The results showed that the six element configuration provided the 
most reduction in run time with a savings of approximately 75% associated with reducing 
the number of elements from 21 to six. The ten and four element models showed a run time 
reduction ranging from 62 to 68%, respectively. 

On the basis of the trade-offs in simulation dynamic characteristics and computer run time, 
a conflguration was selected for further evaluation of run time improvements. Both the bleed 
and ten element configurations showed similar results in terms of transient response time 
and computer run time. Due to the emphasis placed on frequency response improvements, 
the configuration with ten similarly sized elements was selected over the six element (bleed) 
model on the basis of its better frequency’ response characteristics. If model simplification 
or minimum run time had been emphasized in this work, the four or six element models 
would have been better choices. The ten element configuration was used to evaluate further 
run time improvements, using variations in the algorithms for prediction of program iterative 
solutions. 

Predictor Algorithms 

Computation of transients in the SOAPP simulation system requires iterative solution of a 
large number o^ non-linear equations. The program execution time is partially dependent 
on how accurately the solution values of independent variables can be predicted. As the 
accuracy of the predictions improves, the number of iterations required for solution is decre- 
ased, and the computation time reduced. Prediction of the independent variable values from 
solutions at previous time steps is one of the techniques employed in the SOAPP system solu- 
tion routines. 


17 



A study of prediction algorithms was made using the ten element compression system. Four 
different prediction formulas were evaluated u^g step and sinusoidal inlet pressure (F*r2.(P 
disturbances, and a high compressor discharge pressure-time history as input transients. Tlw 
Pt 4 q pressure-time history was derived from a fuel step transient induced into the baseline 
engine simulation (Figure 50). For the purposes of this study, constant values of the time 
increment and the iteration tolerance were employed. The parameter used to evaluate the 
performance of the prediction algorithms was the average number of iterations per solution. 
The four prediction formulas are derived below. 

1 . The simplest prediction scheme is the use of values obtain^ from the most recent 
solution. For the transient case where the effect of the disturbance has died out, 
this predictor offers the best performance since none of the independent variables 
are changing. This is a trivial case but, as will be seen, the response to a step change 
appears to lend itself quite well to this prediction algorithm. 

2. The next level of complexity is represented by a linear predictor of the form: 

X = at+b 

where X = independent variable (9) 


t = time 

If the two past values of ‘X’ & ‘t’ are known, the coefficients ‘a’ and ‘b’ can be calculated 
such that 


/x, -Xj\ /xjt, -X,tj\ 

j' 


( 10 ) 


For the evenly-spaced time interval case (t - 1| - t| - 12), 
This reduces to 


X = 2X| - X 2 


( 11 ) 


3. A quadratic form represents the highest order predictor included in this study 
X = at^ + bt + c (12) 


Again, for equally-spaced time increments, the following formula is obtained 


X = 3Xi - 3X-, + X, 


for t-ti =ti - U = t.,-t 


(13) 


1 M ‘2 ‘2 ^3 

4. An exponential predictor was also investigated in the study. Its form is 


X = ae**^ 


(14) 
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The equally-spaced form of this predictor is 

, ( 15 ) 

X = Xi 2 /X 2 

Table II summaiUea the results of the study. The results show P|t«s ^ 

S1CSwro/the,u»^^^^ 

to be exactly equivalent to the linear. 

This observation is confirmed upon examination of the linearised form of the exponential 
predictor: 

AX / 2 AX,\ AX2\ 


u Kv«Y~Y«Xthis reduces to the linear predictor formula. Most transient 

and the exponential predictors perfoim in a virtually identical manner. 

turbance (^74 q r®*'*^**)- 

stems from the fact that solution p3,^n„o„ luhough the set of 

Sear, whichTso indicates that this type of prediction should perform best 

Based on the results of this study. c°Jes, 

rSh'strpdLti^^ri^:— 

past value predictor. 

TABLE n 


Predictor Algorithm Study Results 



The results obtained from the analysis of the stage-by-stage and reduced order cmnpression 
system simulations led to the following conclusions: 

1 . The frequency response characteristics of dynamic simulations are dictated by the 
largest dynamic elements. The reduced order compression system simulations 
which employed dissimilar element sizing (i. e. fan lumped/compressors staged 
and fan staged/compressors lumped) showed response characteristics similar to 
lower order systems containing fewer dynamic elements of more uniform sizing. 

2. The uniformly sized 10 element compression system model provided the optimum 
configuration for evaluation in the engine simulation because it essentially retained 
the frequency response characteristics of the stage-by-stage configuration while 
offering a significant reduction in computer run time. For other applications 
where model simplicity or minimum computer run time is emphasized, the four 
element or bleed configurations represent optimum models. 

3 ^ The linear predictor algorithm should be used for computation of transient cases, 

because it generally required the smallest number of iterations to reach solutions 
at each time step of a transient. 


Modified Engine Simulation 

The third and final step toward improvement of simulation dynamic characteristcs was 
evaluatipn of the effects of the compression system model modifications on the frequency 
response and transient response characteristics of the engine simulation. To perform the 
evaluation, the baseline engine simulation was modified by incorporating the ten element 
compression system model. The modified engine simulation was used to define gain and 
phase frequency response characteristics and pressure-time histories of response to steps in 
inlet pressure and fuel flow. These dynamic characteristics were compared to results obtained 
from both the baseline engine simulation and NASA engine tests. As a final step in the 
evaluation, the modified engine simulation was subjected to discrete frequency sinusoidal 
inlet pressure disturbances of varying amplitudes to determine the stability limits of the 
simulation. These results were compared to NASA test data for which the inlet pressure 
amplitude at engine surge/stall was defined. 

The following discussion provides a brief description of the dynamic configuration of the 
modified engine simulation. The dynamic characteristics of the simulation are presented and 
compared to results obtained from the baseline engine simulation and NASA engine tests. 
Tlie dynamic characteristics are discussed in terms of frequency r'^sponse (gain and phase), 
transient response (pressure versus time), and stability limits. 

Configuration 

The baseline engine simulation was modified by implementing the ten element compression 
system model and linear predictor algorithm selected from the reduced order analysis. The 
resulting dynamic configuration of the modified engine simulation is shown in Figure 5 1 
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compared to an engine side-view. The number of dynamic elements was 27. compared to 
20 in the baseUne simulation. This increase came from six addition^ elements m the fan, 
low and high compressors plus separation of the high compressor discharge dafuser from 
the compression system. The state variable order of the system was 83, which included 
mass, momentum, and energy for each of the 27 dpamic elements plus fuel flow and core 
exhaust area. The simulation is described in detail in Appendix A. 


Dynamic Characteristics 

The dynamic characteristics of the modified engine simulation were established in three 
formats. First, the frequency response at stations 2.3, 2.3F, 3.0 and 4.0 was demed using 
both the transfer function and step response methods. The results of these metoods were 
compared to both engine test and baseline simulation results. Second, time lustones of 
response to steps in inlet pressure and fuel flow were computed for stations 2.3, 2.3h, 

3.0 and 4.0 and compared to baseline simulation results. Third, discrete frequency 
pressure sine wave transients of varying ampUtudes were induced to determine the stability 
limits of the simulation. Results were compared to engine test data where the pressure 
amplitude at surge/stall was defined. 

Frequency Response 

The main objective of this program was to develop an engine simulation which showed im- 
proved frequency response characteristics to at least 100 Hz. To accomplish this, the fre- 
quency response characteristics of the simulation were defined using both the step response 
and transfer function methods. The results obtained from the modified engine simulation 
were compared to engine test results and baseline simulation results to assess improvements 
in simulation frequency response characteristics. In addition, frequency response character- 
istics derived by the two methods were compared to each other to determine which gave 
the better representation of the actual engine. 

The frequency response characteristics computed for the modified engine simulation are 
shown in Figures 52 through 67 which show response computed by the transfer function 
and step response methods compared to the engine test data and transfer function 
characteristics of the baseline engine simulation. Calculation of transfer functions for the 
modified engine simulation was found to be more difficult than was experienced with the 
baseUne engine simulation. Implementation of variable scaling of the program iteration 
tolerances and linearization of performance characteristics were required 
transfer function calculations. With the variable scaling, all iteratnns were stdl not withi 
desired tolerances, however, the frequency response results appeared to provide sufficient 
accuracy for some comparisons. The effects of solutions outside of tolerance bands were 
most evident in the core stream response characteristics at condition OPl and in the bypass 
stream characteristics at both OPl and OP2. The most common evidence was unusually low 
gain levels at frequencies up to 10 Hz or a 360® shift in the phase angle. 

The influence of the modified compression system model on the engine simulation core 
stream gain characteristics is shown in Figures 52 through 57 which show that the simulation 
gain levels were affected at frequencies above 50 Hz. The transfer function and step 
response gains showed a continuation of roll-off above 50 Hz which was not present in the 
baseUne engine simulation. At 200 Hz the modified simulation gain levels at station 4.0, 

3.0, and 2.3 were reduced by an average 1 8, 9, and 4 decibels respectively. 



The transfer function gain levels at condition OP2 (Figures 52 through 54) showed im- 
proved resolution of the simulation high frequency response characteristics. The gain re- 
sponse showed resemblence to the NASA test results but appear to be displaced to higher 
frequencies. The resemblence was characterized by a gain peak which occurred at 80-90 
Hz in the simulation, compared to 30-35 Hz in the test data. The simulation gain levels 
predicted by step response at condition OP2 did not resolve high frequency dynamics as 
well as the transfer function. 

At condition OPl, core stream transfer function results showed the effects of program 
solutions reached outside the desired iteration tolerance. The gain levels defined at the 
high compressor discharge (Figure 55) disagreed with step response results over the range 
of 1 to 40 Hz. Agreement was good over this frequency range at condition OP2. At 100 Hz, 
the transfer function predicted a gain peak at station 4.0, but no similar peak was seen 
in the gain response at stations 3.0 and 2.3 (Figures 56 and 57). These results were 
inconsistent with condition OP2, NASA test data, and fan bypass gain response in the 
baseline engine simulation. At condition OP2, simulation gain peaks at the three core 
stations were aligned at a frequency of 80-90 Hz. At both OPl and OP2 the NASA test data 
showed alignment of peaks at 30 to 35 Hz and 80 to 94 Hz. All of these examples showed 
influence of downstream components on the gain characteristics. The gain levels at stations 
2.3 and 3.0 were showing influence of the burner and high compressor dynamics. The fan 
bypass gain peaks shown by the baseline simulation indicated that response was influenced 
by the fan duct dynamics. 

The step response gain levels at condition OPl (Figures 55 through 57) indicated the same 
results as condition OP2. The gain characteristic rolled-off to 200 Hz with no consistent pre- 
diction of gain peaks. The influence of the compression system modifications on simulation 
core stream phase is shown in Figures 58 through 63. These curves show that the phase angles 
were affected over the full spectrum of frequency ( 1 to 200 Hz). In general, the transfer func- 
tion predicted more lag for the modifled engine simulation than was predicted for the base- 
line simulation up to the frequency where phase peaks were calculated. 

At condition OP2 (Figures 58 through 60) the transfer function phase showed improved 
resolution of the high frequency response characteristics of the modified simulation. Peaks 
in phase were consistently predicted at all three core stream stations in a frequency range 
of 70 to 76 Hz. The phase characteristics showed some resemblence to NASA engine data 
although again appear to be displaced to higher frequency. The NASA data showed phase 
peaks at approximately 35 Hz. The NASA test data also showed less lag at the phase peaks. 

The core stream phase computed at condition OP2 by the step response method (Figures 58 
through 60) showed less lag than transfer function phase; this was consistent with results 
from the baseline simulation. As in the case of the modified simulation gain, there was no 
consistent prediction of phase peaks with the step response method. 

The core stream phase at condition OPl is shown in Figures 61 through 63. The high 
compressor discharge station (Figure 61) again showed the effect of reaching program 
iterative solutions outside of desired tolerances. The effect was indicated by comparison of 
transfer function and step response results. Over the range of 1 to 75 Hz, the transfer 
function phase showed less lag than step response which was counter to othei simulation 
results. The transfer function phase at the high compressor discharge alio showed a steep 
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roll-off in phase angle above 65 Hz which gave some resemblence to engine data at 
approximately 25 Hz. No corresponding phase response was shown at the upstream stations 
(3.0 and 2.3), again indicating that the influence of dynamic characteristics of downstream 
components was not being resolved at the upstream stations. 

The step response phase computed for the modified engine simulation at condition OPl 
(Figures 61 through 63) indicated results similar to condition OP2. With exception of the 
high compressor discharge station, the step response phase showed less lag than transfer func- 
tion phase and made no consistent prediction of phase peaks or steep roU-off. 

The frequency response characteristics of the fan bypass stream discharge (station 2.3F) 
are shown in Figures 64 through 67. Comparison of the transfer function gain levels from 
the baseline and modifled simulations (Figures 64 and 65) showed that response peaks 
continued to be predicted by the modifled simulation, although the levels were generally 
lower than the baseline simulation above 50 Hz. (The -20 db level at 9 Hz predicted by the 
modifled simulation transfer function at condition OP2, (Figure 64) was probably a result of 
calculation inaccuracy, not system characte.i-tics). The gain levels computed from step re- 
sponse generally showed the same results as the baseline simulation; response peaks above 30 
to 35 Hz were significantly attenuated. 

Comparison of transfer function phase results in Figures 66 and 67 show that the modified 
simulation phase angles showed less lag than the baseline and were closer to engine test levels. 
(The 360 degree phase shift shown at OPl was again due to calculation in accuracy). The step 
response phase indicated less lag than the transfer function predictions, consistent with vir- 
tually all other results. 

Transient Response 

Pressure-time histories at stations 2.3 F, 2.3, 3.0 and 4.0 are shown in Figures 68 through 75 
for a 1.0% step in inlet pressure. These figures compare results obtained from both the baseline 
and modified engine simulations. Figure 68 shows that the station 2.3F fan bypass stream 
discharge pressure at condition OPl exhibited a stairstep characteristic of response resulting 
from back and forth disturbance propagation in the fan ducts. Results were similar for both 
the baseline and modified simulations. 

Pressure-time response in the core stream at condition OPl shown in Figures 69 through 71 
reveals the response of the baseline and modified simulations was different. The baseline 
simulation core stream stations responded in stairstep fashion similar to the fan bypass 
stream while the modified simulation did not. These results indicated that the baseline 
simulation coupled the fan bypass and core stream dynamics while the modified simulation 
did not. Similar results were shown at condition OP2 as illustrated in Figures 72 through 75. 

The pressure-time response to a 1 .0% fuel flow step is shown in Figures 76 through 79. 

Figure 76 compares at condition OPl, the baseline and modified simulation response at the 
high compressor and fan bypass stream discharge stations (4.0 and 2.3 F). The results at these 
stations were similar for both simulations where the pressure at station 4.0 rose sharply over 
the first 0.01 seconds while the pressure at 2.3F showed significant rise at about 0.02 seconds. 


23 



The delay in station 2.3F response indicated propagation of the pressure wave resulting from 
the fuel pulse rearward into the exhaust system and then forward through the fan bypass 
ducting to the fan bypass stream discharge. 

The pressure-tim« response at condition OPl is shown for core stream stations 2.3 and 3.0 
in Figure 77. At these stations, results from the two simulations were different. The baseline 
simulation showed a small rise in pressure over the first 0.01 seconds followed by a second 
larger rise in pressure starting at approximately 0.02 seconds. These characteristics indicated 
forward propagation of the pressure wave resulting from the fuel pulse over the first 0.01 
seconds followed by coupled response to the fan bypass beginning at approximately 0.02 
seconds. The modified simulation response showed tnat the fan bypass-core coupling was 
not present; only forward propagating response to the fuel pulse was seen over the first 
0.01 seconds. The results at condition OP2 were the same as shown in Figures 78 and 79. 

The differences in fan coupling characteristics of the baseline and modified simulations were 
attributable to the respective simulation methods. The fan representation used in the 
baseline simulation, provided coupling of the bypass and core stream pressures through the 
performance characteristic illustrated in Figures 6 and 7. The modified simulation contained 
no coupling. The performance characteristics of the fan bypass and core stream were 
described in terms of separate 0, 0 and X characteristics with the splitting streamline 
computed as a function of bypass ratio. In developing the modified compression system 
simulation, attempts were made to compute the splitting streamline during transients, 
however, the use of a constant flow per unit area assumption in the computation appeared 
to give dynamically unstable results. Transfer functions computed for the simulation showed 
positive roots of the characteristic equation and confirmed tne dynamic instability of the stream- 
line calculation. As a remedy, the streamline location was computed at the initialization 
point of the transient and held constant for the duration of the transient. As shown by the 
pressure and fuel flow step transient results, coupling therefore was not present in the modi- 
fied simulation. 

Stability Limits 

The stability limits were assessed as a final step in evaluation of the dynamic characteristics 
of the modified engine simulation. At condition OPl, inlet pressure sign waves of varying 
amplitude were induced into the simulation to determine the approximate amplitude at 
which simulation computations became unstable. The results were compared to engine test 
data where the sine wave amplitude at engine surge/stall was defined. 

A comparison of simulation and test results (Figure 80) show that the simulation instabilities 
occurred at wave amplitudes at least twice as large as the amplitudes which induced surge/ 
stall during the engine test (60% to 70% for the simulation versus 30% for the engine). The 
instability induced in the simulation occurred in the fan bypass stream. This result is 
illustrated in Figure 81 where the transient path of operation is shown on the fan map for 
both stable and unstable cases at 10 Hz. The stable case shows orbiting of the transient 
operating point in a clockwise direction, compared to the unstable case which shows 
a rise in pressure followed by a large reduction in airflow to the point the simulation 
computation became unstable. The airflow characteristic shown for the unstable case is 
typical of test results experienced by Pratt & Whitney during engine surge cycles. 
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excuisions (Figure 83) were weU below the steady state surge Une. 
TheresuiUof^eNASA^e^-^^^^ 

in - 0^0—- - rorm o( 

routing stall, and culminated with high compressor surge. (See Reference 5). 

Tire comparison of simulation and engine stabUity limits led 
me simition predicted the fan to be '7' “!*" 

Sud^"dS“Suln 

required to properly predict engine core stream dynamics. 

The resulte which were obtained from frequency response, time history response, and stabUlty 
oT^sLulation suggest areas where further work is required to unpnwe the 
staulaUon dynamic characteristics. These areas are identified under the discussion of results, 
conclusions, and recommendations sections of this report. 
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DISCUSSION OF RESULTS 

The investigations made during this program produced the following key results. 

I . Use of the dynamic energy equation in combination with dynamic mass and 
momentum equations is required to properly simulate the propagation of 
aerodynamic disturbances through the system. With the energy equation 

^.delayed response of downstream pr^ures and ’ 
emperatures with inlet induced disturbances; without the energy equation 
downstream response was immediate and therefore improper. ^ 

^ thTt Tm rin compression system simulation showed 

^u d ’ characteristics of a stage-by-stage simulation 

^u d be obtained with a more simplified dynamic conHguration while 
realizing substantial savings in computer execution time. The analysis showed 
S elements should be uniformly sized to achieve the best 

•°'"cr frequency limitations of iJ^e elements 

hat compression systems never require simulation on a stage-by «tage basis 

s^eeSroWem of interest in a ' 

specific problem. As the maximum frequency of interest becomes higher the 

required number of dynamic elements should increase in a manner wS 

aintains a sim^ilar relationship between element size and propagation time at 

the maximum frequency. The real significance of the reduL order andys!s 

was^e approach that was used. By example, it illustrated a method which 

«dth*TesMcr? 7*^*’**"^ ^ simulation dynamic configuration which is optimum 
with respect to frequency response and program execution time. 

algorithm analysis showed the linear form to produce the best 
erall efficiency in computation of simulation transients. The good performance 
of the algorithm resulted from its characteristics of minimal amplification of 
errors associated with program iteration tolerances and the tendency toward 
hneanty associated with the small time steps used in transient computations. 

^e transfer function provided better resolution of simulation high frequency- 
response characteristics than the transient methods. This was characterized by 
the peaks ga.„ and phase frequency response curves. In cases whe4 tte 
ransfer functions were calculated with satisfactory accuracy, core stream 

consistently resolved at all stations and aligned at constant 

\ ^ ^ {unction appeared to resolve 

sponse peaks associated with the dynamics of the fan downstream ducting. 

The frequency response characteristics obtained from the transfer function and 
^ulation transients may differ because of the nature of their respective solutions 
I ransfer functions are computed by direct solution of the linear system equations 
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derived from the nonlinear simulation. Errors result in this method from the itera- 
tion tolerances needed to compute the linear system model and from the numerical 
evaluation of the transfer function which requires raising of numbers to large powers. 
Transient computations are not free from error because of the tolerances required for 
iterative solutions at each time step of the transient. The tolerances introduce pseudo 
dynamics into the transient which are reflected in the frequency response characteris- 
tics. 

The differences inherent in the two methods of defining simulation frequency 
response characteristics suggest that the transfer function method can be used to 
advantage. For example, transfer functions can be used to measure how well 
simulation transients reflect system response characteristics. The advantage of 
the transfer function also endorses current applications of the method for 
analysis of control system operation. 

5. Modifications to the compression system model primarily affected frequency 
response characteristics of the engine simulation above 50 Hz. The effects 
were most clearly illustrated by the reduced gain levels which showed continued 
roll-off above 50 Hz. The transfer function response showed gain and phase 
peaks similar to test data, indicating that the dominant dynamic characteristics 
of the engine were simulated but appears to be displaced to higher frequencies. 

The low frequency response characteristics of the engine simulation were essen- 
tially unaffected by the modifications to the compression system model. In 
an absolute sense, close agreement between the modified simulation and test data 
remained confined to a maximum frequency of 10 to 20 Hz in the core stream 
and 40 to 65 Hz in the bypass stream. 

The simulation frequency response characteristics showed that, although modifi- 
cations to the compression system model improved resolution of the dominant 
dynamic characteristics of engine, the approach used in this program was not 
sufficient to achieve close agreement with test data over the frequency spectrum 
of 0-200 Hz. The approach used in the program was centered around the 
influence of compression system dynamic element sizing and results showed that 
further modifications to the engine simulation are required. Three possibilities 
can be suggested at this point. 

The reduced order compression system analysis showed that dynamic element 
sizing should be uniform to achieve the best frequency response characteristics. 
When the ten element compression system model was incorporated into the engine 
simulation, the dynamic element sizing of other components was not reduced 
for compatibility with the compression system. This would suggest that if further 
modifications were pursued, dynamic element sizing of other engine components 
should be reduced and evaluated. 

The comparisons of simulation and test frequency response diaracteristics at the 
fan bypass station showed agreement at higher frequencies (40 to 65 Hz) than at 
core stream stations (10 to 20 Hz). The bypass flow system component compo- 
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sition is dominated by ducting while the core stream is more heavily composed 
of turbomachinery. These facts could suggest that the dynamics of the ducting 
were more satisfactorily modeled than the turbomachinery dynamics. The simu- 
lations used in this program placed the boundaries of the compression system dy- 
namic elements in the gaps between rotor and stator rows and assumed linear dis- 
tributions of the flow properties between the element inlet and exit boundaries. 

This assumption therefore did not account for the reduction in flow area and 
associated increase in Mach number experienced by the flow as it progresses 
through the airfoil rows. The effect of the assumption used in this investigation is 
to compute Mach numbers through the compressors which on the average are too 
low and in turn make propagation time of disturbances traveling upstream too fast 
This would suggest that the description of flow property distributions through the 
turbomachinery should be modified and the effects on frequency response evaluated. 
The engine simulations assembled in this program were composed of dynamic 
elements representing only the engine components. Testing of the actual engine 
was performed in an altitude test facility which included inlet ducting of signifi- 
cant length ahead of the engine inlet plane. This ducting potentially influenced 
the dynamic characteristics of the engine-test facility system. It is expected 
that the inlet ducting would shift the resonances of the system to lower frequen- 
cies, therefore if the ducting ahead of the engine was included in the simulation, 
the agreement between simulation and test frequency response characteristics 
might be improved. Addition of the inlet ducting to the simulation is a modifi- 
cation which should also be considered. 

6. Comparison of the step transient results from the baseline and modified engine 
simulations showed that coupling of the fan bypass and core dynamics have a 
significant effect on core stream transient operation. These effects were illustra- 
ted by comparison of the results from the baseline and modified engine simula- 
tions. Coupling was present in the baseline simulation where transient results 
showed that the core pressures responded to disturbances propagated forward 
through the fan bypass ducts. No coupling was present in the modified simula- 
tion, and core response to fan bypass stream disturbances were not shown by 
simulation transient results. 

These results are significant to the engine design process and underline the 
importance of the role of the dynamic engine simulation. The results show that 
dynamic characteristics of the turbofan bypass flow system can potentially play 
a significant role in the transient operation of the core stream compression 
system and thereby influence the stability of operation. Design characteristics 
such as length of the fan ducts, proximity of the flow splitter to the fan, and 
use and placement of fan airfoil shrouds may all play an important role in 
setting the compression system stability limits. The capability to accurately 
simulate the effects of these design characteristics is therefore valuable to the 
design process. 

7. The stability limits of the modified engine simulation were more than twice that 
shown by engine test results. The simulation results showed that instability 
originated in the fan bypass stream while test data showed origin in the core 
stream. 
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The comparison of simulation and test stability limits also showed that further 
refinements to the simulation are needed to achieve a icUable estimate of the 
engine stability limits. Further revisions to 

dvMmic element property distribution were proposed above. Fw byp^/ro 
S^uS waTalso proposed as potentiaUy very important to stabthty of ojrera- 
«r S cSd«aUon to be mad. is that instabUities usually o"!”"'.'" 

nwdeling used during this program might need expansion to account for the cir 
cumferential properties of rotating stall. 

A variety of tests was required to reveal important dynamic 

the simulation. The dynamic characteristics of the simulations assemble u g 

tt^rTwere evaUted using three different tests: 1) frequency respons^ 

S^n a^d phase versus frequency). 2) transient response pressure • me). 

aSd 3) stabdity limits (Pt 2.0 amplitude). Results of these tests illu 

trated by example the importance of each test. 

The frequency response characteristics provided an overall view of the Emulation 

dynamics across the frequency spectrum, and suggested 

further refinements to the engine simulabon in areas 

sizing and distribution of properties across the f 

frequency range of 0-20 Hz, where the simulation stability was evaluated, 

frequency response characteristics of the baseline engine simulation, the 

modified engine simulation, and test data where in closest agreement. 

The transient response characteristics revealed fan bypa^co^^ 
coupling was present in the baseline engine simulation but * 

modified engine simulation. The results illustrated that the coupling effect 
has a potentially important effect on engine stability. These results were not 
apparLt from the frequency response characteristics of the two simulations. 

The stability check provided the final test of simulation accuracy and showed 
m«"simulation termemenU are necessary to adtieve a rehablc predrebon 

of engine stability limits. 

In the process of developing the simulations used in this program, three problems 
of a numerical nature were encountered. Two of the problems were assoaated 
with computation of frequency response by the matrix method; the third problem 
was associated with stability of the model formulation. 

The fust attempt at computing frequency response characteristics by the matrix 
method gave results which showed non zero gain levels at frequencies near zero 
Hz. These results were derived from finite difference perturbations which were 
one sided; i.e.. all in a positive or negative direction. Experiments with the per- 
turbaUon size and iteration tolerances showed that the low frequency offsets in 
gain level persUted while results at higher frequencies were repeatable. Further 
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tests showed that by changing the direction of the perturbation, i.e.., from posi- 
tive to negative, the direction of the low frequency gain offset also changed 
while higher frequency results again were repeatable. On the basis of these tests, 
the fmite differences routines of the program were modified to compute the 
matrix elements from perturbations in both the positive and negative directions. 
This modiHcation virtually eliminated all offsets in low frequency gain levels as 
long as all program iterations were completed within input tolerance levels. 

All of the frequency response characteristics derived by the matrix method for 
the baseline engine simulation were successfully computed within the iteration 
tolerance bands. Two cases computed for the reduced order compression system 
models reached solutions with some iterations falling outside the prescribed 
tolerances. Further difHculty was encountered in computing matrices for the mod- 
ified engine simulation. Scaling of the program tolerances was required to reach 
completed solutions, and then some iterations were still outside of tolerances. 

The effects of reaching solutions outside of prescribed tolerance bands were 
usually shown by offsets in gain and phase levels in the low frequency range. In 
addition, the case run time was also extended, reflecting the attempts of the 
program to reach satisfactory solutions. The causes of these iteration problems 
could not be determined during this contract, however it was observed that the 
occurences correlated with simulation configurations containing dynamic elements 
of dissimilar sizing. 

The final numerical problem encountered with the simulations was associated 
with the dynamic formulation of the splitting streamline in the improved com- 
pression system model. An assumption of constant airflow per unit area was 
used to compute the radial location of the streamline separating the core and by- 
pass flows of the fan. Initial attempts to compute streamline locations during 
transients appeared to be unstable. In this case, the matrix computations proved 
to be very helpful to the problem because the characteristic equation showed 
positive roots and conflrmed instability. To stabilize transient computations, 
the simulation was modified to compute streamline location at the initialization 
point and held constant for the duration of the transient. 



CONCLUSIONS 


The results obtained from investigation of the dynamic characteristics of the mixed flow 
turfoofan engine digitally simulated under this program have led to the following conclusions: 

1. High frequency simulations should incorporate dynamic forms of the mass, 
momentum and energy equations. 

2. High frequency engine simulations using finite dynamic elements can be assembled 
in a manner where the dynamic configuration is optimum with respect to dynamic 
characteristics and computer execution time. The dynamic configuration 

should employ dynamic elements of uniform length. The number of required 
elements will be dictated by the maximum frequency of concern. 

3. Implementation of a linear predictor algorithm improved the execution time of 
the simulations investigated under this program. 

4. The transfer function method provides good definition of system frequency re- 
sponse characteristics. The transfer function method can therefore be used to 
assess the frequency response characteristics derived from simulation transients 
and, as a corollary, should not be interpreted as necessarily representative of the 
frequency response characteristics of simulation transients. 

5. Dynamic coupling of the fan bypass and core streams can have a significant in- 
fluence on core stream operation and stability. Accurate prediction of these 
effects can be important to the engine design process. 

6. Tests of the simulation frequency response, transient response, and stability are 
all useful in assessing the simulation characteristics. The frequency response char- 
acteristics provide an overall view of the dynamic characteristics of the simulation 
and define the valid frequency range when compared to actual data. The transients 
are useful in identifying the dynamic coupling between the system components. 
The stability test provides a final check of the simulation when compared to test 
data, establishes the credibility of predicted engine stability limits. 

7. Resizing the compression system finite elements improved the dynamic character- 
istics of the engine simulation above SO Hz but showed that additional refinements 
are required to obtain close agreement between simulation and actual engine dyna- 
mic characteristics. The modified engine simulation showed good bypass stream 
dynamics to 60 Hz and good core stream dynamics to 20 Hz. 
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RECOMMENDATIONS 


Further development of the turbofan engine simulation assembled during this program 
is required to obtain close agreement with actual engine dynamic characteristics. On the 
basis of the findings in this program, it is recommended that refinements be made in the 
following areas: 

1. Formulation and evaluation of equations which properly model the dynamic 
coupling of the fan bypass and core streams. With these revisions to the fan 
dynamics characteristics, it is possible that another configuration of the reduced 
order compression system would be optimized with respect to frequency response 
and computer run time. 

2. Resizing of the dynamic elements of other engine components for compatibility 
with the reduced order compression system model. The frequency response char- 
acteristics of the simulation would be improved but at the expense of increasing 
the number of dynamic equations. 

3. Redefinition of the flow properties across the dynamic elements to recognize the 
area variations which occur as flow passes through the gaps and airfoils of the 
turbomachinery. These modifications would probably introduce additional com- 
plexity to the simulation. 

4. Model the inlet ducting to assess its effect on frequency response. This would 
also increase the number of dynamic equations. 

5. Modeling of rotating stall may be needed but should be implemented as a last 
resort. Test experience with the TF30 engine has shown low compressor rotating 
stall to be a common form of initial instability in the compression system. This 
implies that modeling of this process could finally be required to obtain agreement 
between simulation and test stability limits, however, this approach should be a 
last resort because of the substantial complexity that would be added to the simu- 
lation. 
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Figure 2 Configuration of Baseline Engine Simulation 
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Figure 6 Fan Core Pressure Ratio Map 


Figure 7 Fan Core Temperature Ratio Map 
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Figure 8 Effect of Energy Equation on Station 7. OF Pressure and Temperature Response 
to a Step in Inlet Pressure 
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Figure 9 Illustration of Wave Artalysis for Evaluation of Gain and Phase 
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Figure 12 Fan Core Stream Cain at Station 2.3 
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Figure 13 Fan Core Stream Phase at Station 2.3 
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Figure 15 Low Compressor Phase at Station 3.0 
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Figure 16 Fan Bypass Stream Gain at Station 2.3 F 


lOOr 


CONDITION OP1 



€0 

m 



- 150 - 

- 200 - 


NAS* test data 

. 1 -^ — Baseline engine simulatioi; transfer function response 
Baseline engine simulation step response 
# Baseline engine simulation discrete frequency response 


fl 1 1 1 1 1 1 i 1 1 1 

1 1 1 1 1 1 1 11 

h^bbuiiii^^b^bi 


0.1 1 

10 

100 



Frequency ~ Hz 


Figure 17 Fan Bypass Stream Phase at Station 2.3F 
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Figure 1 9 High Compressor Phase at Station 4. 0 
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Figure 20 Fan Core Stream Gain at Station 2.3 
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Figure 21 Fan Core Stream Phase at Statio't 2.3 
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Figure 25 Fan Bypass Stream Phase at Station 2.3F 
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Figure 26 Compression System Simulation, Stage-By-Stage Configuration 
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Figure 27 Bypass Ratio Effects on Dividing Streamline Location and Core Stream Stator 
Operation 
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Figure 28 Airfoil Lift Characteristics 



Stall 

limit 


Figure 29 



Iniet air angie 

Compressor Airfoil Loss Characteristics 



Lfgend: 

“ Step response 
Transfer function response 


I lie 


1 


ON 

o 


CONDITION 0P1 







4ir 




'«J» v*-,;o. 

































s 


Fan 



0 20 40 60 80 100 120 CHI 



Figure 35 Reduced Order Cortipression System Simulation, Ten Element Configuration 
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Figure 38 Reduced Order Compression System Simulation, Fan StagedlCompi 
Lumped Configuration ^ 
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Figure 39 Compression System Simulation Gain at Station 4.0 
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Figure 40 Compression System Simulation Phase at Station 4.0 
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Figure 42 Compression System Simulation Phase at Station 4.0 
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Ffgurt 43 Correlation of Compression System Gain With Number of Dynamic Elements 
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Figure 44 Correlation of Compression System Phase With Number of Dynamic Elements 
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Figure 45 High Compressor Pressure-Time History of Response to a 1.0% Step in Inlet 
Pressure 
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Figure 46 High Compressor Pressure - Time History of Response to a 1.0% Step in Inlet 
Pressure 
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Figure 49 Compression System Simulation Computer Run Time for Transfer Function 
Calculation 
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Figure 51 Configuration of Modified Engine Simulation 
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Figure 52 High Compressor Gain at Station 4.0 
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Figure 54 Fan Cote Stream Gain at Station 2.3 



Pt4.o gsin ^ decibels 


CONDITION OP1 


NASA test data \ 

Modified engine simulation transfer function response 
Modified engine simulation step response 

Baseline engine simulation transfer function response 


Frequency 


Figure 55 High Compressor Gain at Station 4.0 
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Figure 56 Low Compressor Gain at Station 3.0 
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Figure 59 Low Compressor Phase at Station 3.0 
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Figure 60 Fan Core Stream Phase at Station 2.3 
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Figure 62 Low Compressor Phase at Station 3.0 
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Figure 63 Fan Core Stream Phase at Station Z 3 
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Figure 65 Fan Bypass Stream Gain at Station 2.3F 
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Figure 66 Fan Bypass Stream Phase at Station 2.3F 
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Figure 67 Fan Bypass Stream Phase at 
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Figure 69 Fan Core Stream Pressure - Time History of Response to a 1.0% Inlet Pressure 
Step 
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Figure 71 High Compressor Pressure - Time History of Response to a 1.0% Inlet Pressure 
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Figure 74 Low Compressor Pressure - Time History of Response to a 1.0% Inlet Pressure 
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Figurs. 75 High Compressor Pressure - Time History of Response to a i.0% Inlet Pressure 
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Figure 77 Fan Core Stream and Low Compressor Pressure • Time Histories of Response to 
a i.0% Fuel Flow Step 
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APPENDIX A 

DESCRIPTION OF MODIFIED ENGINE SIMULATION 


INTRODUCTION 

This appendix describes the modified dynamic engine simulation evaluated during this pro- 
gram. This appendix is divided into three sections which are 1) performance characteristic 
used to simulate each component, 2) dynamic equations employed, and 3) a discussion of 
the program solution techniques. 

The method of defining performance characteristics of each engine component is described 
in the Component Performance Characteristics section. The components described are the 
compressor, burner, turbine, ducts, exhaust nozzle and the control. The mdividual compo- 
nent characteristics describe the steady state operation of each component. 

The Dynamic Equation section describes the method of simulating the dynamic character- 
istics of each component. The method uses dynamic forms of the continuity, momentum, 
and energy equations. The form of the equations used in the simulation is presented. Two 
dynamic computations, rotor dynamics and dynamic mixing, do not directly use the three 
basic equations to describe engine operation. The method of implementing rotor dynamics 
and dynamic mixing is also presented. The use of the dynamic equations in conjunction 
with the steady state performance characteristics is described. 

The Solution Techniques section discusses the methods used to solve the simulation equa- 
tions. The form of the simulation equations is implicit and iterative (relaxation) methods 
are used to arrive at solutions. 


COMPONENT PERFORMANCE CHARACTERISTICS 


Compressor 

Description 


The TF30-P-3 compression system consists of a three stage fan, a six stage low pressure 
compressor and a seven stage high pressure compressor. Figure A-1 presents a cross section 
of the TF30-P-3 compression system where key stages have been identified. The fan com- 
presses all of the air coming into the engine. Part of the air flows through the core, i.e., 
into the low and high compressors, to be compressed to a higher pressure. The remaining 
flow goes only through the fan and is bypassed around the core. The ratio of the bypass 
flow to the core flow is the bypass ratio. 

The compression system used in the modified engine simulation was based on individual 
row performance characteristics. A row consisted of either a rotor or a stator and a gap 
behind the rotor or stator. The rows were interconnected to form the entire compression 
system. 

The compression system pressure and temperature change were defined through the use of 
characteristic maps for each individual row. The temperature rise was calculated from total 
temperature maps and the pressure change from static pressure maps. Total and static 
pressures in each row were related by standard compressible flow relationships. 

Mathematical Model 


In a row-by-row model, parameters arc computed at the entrance to each stator or rotor 
and at the entrance to the gap (exit of the rotor or stator). The exit performance parameters 
of each row becomes the entrance conditions of the following icw. Performance modifiers 
are used at various places in the compression system to account for off design operation of 
the system. These include: 

1 . Modification of the fan bypass stream and fan core stream performance charac- 
teristics to account for bypass ratio changes. The characteristic maps for the fan 
are determined for a design bypass ratio (BPR). The modifier allows calculation 
of fan performance parameters at bypass ratios other than design. 

2. Modification of the performance characteristics in the vicinity of the splitter to 
account for flow separation at the splitter and the resulting change in flow angle. 

Inlet Guide Vane 


The normal simulation input parameters at the engine entrance are total pressure, total 
temperature and air flow. The row performance characteristics are based on static pressure 
therefore a static pressure must be calculated prior to starting the first row (inlet guide 
vane) computation. To determine the static pressure from the given parameters, a flow 
parameter is calculated as 


F^ENTRANCE * 


^a ENTRANCE >/ ^T ENTRANCE 
**T ENTRANCE ^ENTRANCE 


(Al) 
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The static to total pressure ratio is determined from the compressible flow relationship 
dtown in Figure A-2. 

The static pressure is calculated as 


**T ENTRANCE 

**S ENTRANCE = /T \ ^ ^ 

(*^T ENTRANCE ™S ENTRANCE ) 

Row Computations 

The entrance conditions to each row are obtained from the exit of the previous row or 
from the inlet guide vane entrance calculations. The rows are interconnected, or stacked, to 
provide a simulation of the complete compression system. The computations performed in 
each component of the row are diagrammed in Figure A-3. The remainder of this section is a 
presentation of the formulas which are used in the calculation of the thermodynamic pro- 
perties through the compressor row. 


The compressor row characteristics are presented as pressure and temperature functions 
Mdiich are correlated with a flow function and the row inlet corrected speed. The flow 
function for a given row is 


«n 


^ag n-1 V^Tg n-l^STAND 
^**Tgn-l^STANDjA„ 


/ \ 

DESn 

n-i/TsTANDJ 

n-l^STANDj 

n 


MOD 1 


MOD 2 
(A3) 


Where 


^N/v/Tjg n-l/^STAND^ 


DES n 


(A4) 


is the design speed for each row and is a constant value determined when the row character- 
istics are generated from test data. 

MODI in equation (A3) is a modification applied to 


^AATg n-l^STANp) DES n 
NVTjgn-l/TsTANDjn 


(A5) 


IIS 


of the fan rows from station 2.0 to station 2.3 and 2.3F shown in Figure A-1. The modifica- 
tion is an equation which is a function of the geometric fan dimensions and the bypass 
ratio. The modiHcation for the core stream is 


MODI » 


(0D„2 + 3ID„2)(1 +bPR) 
2(0D„2 + (1 +2BPR) ID„2) 


1/2 


(A6) 


The modification for the bypass stream is 

(3 00^2 + id„2) ( 1 + bPR) * /2 

2 ^2 + BPR)ODn2 + iD„2^ 


MODI H- 


(A7) 


This modification is made to allow use of the row characteristics at a bypass ratio other 
than the bypass ratio at which the characteristics were generated. The characteristics are 
generated for a given position of the streamline separating the bypass flow and the core flow 
which is a function of the bypass ratio. However, when the bypass ratio changes, the stream- 
line between the bypass and the core stream moves which in turn changes the effective rotor 
or stator length. The change in effective rotor ot stator length in turn eflects the average 
wheel speed. The flow velocities through the rows are proportional to wheel speed and 
MODI can therefore be applied to the speed ratio term of the flow function to the account 
fcM' the effective length variations. 


MOD2 is a modification of the area at core stator 3 shown in Figure A-1. The modification 
is applied at low bypass ratios to account for a reduction in the flow area when flow separa- 
tion exists at the splitter. The modification is 


M0D2 = 0.38462 BPR + 0.5 1 538 


(A8) 


and is limited to a maximum value of 1.0. 

The pressure function {^) and the temperature function (X) are determined from the row 
characteristic maps shown schematically in Figures A-4 and A-5, respectively. 

The stator exit total temperature Ls assumed to be equal to the stator entrance total 
temperature because no work is input to the flow through the stators. 

The rotor exit total temperature is computed using the relationship 
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^Tcn = T'Tgn-l 


1.0 + 


MOD 1 


(NVTjg n-l/TsTANlOn 

2 

(N/s/Tjg n-l/^STANDfDES n 



(A9) 


The rotor and stator exit static pressure is computed using the relationship 


^Scn ^Sgn-1 


1.0 + 


^^ + MOD3 


MOD 1 


^AATgn-l/TsTAND)n 


^V^Tg ii-I^STANd) DES n 


(AlO) 


where MODI is the same value and applied at the same location as the calculation. 
MOD3 is applied to the bypass stream rotor 3, core stream rotor 3 and core stream stator 
3. These locations are identified on Figure A-1 . This modification is applied to account 
for the radial bending of the streamlines in the vicinity of the splitter. 


MOD3 is a dimensionless parameter defined as 


MOD3 = 


MODI 


*^M * “n zg n-1 
- l O-M^zgn-l - 


■( 

^VT'Tgn-l^STAND)DESn ] 

A 


In J 

|NVTTgn-l/^STAND; 


2 

(All) 


where Kj^ is defined in a different manner for each of the three locations where MOD3 is 
applied. 

For bypass stream rotor 3 


^83 BYPASS/^S3 CX)RE 

1.0 + BPR|^3 

(A12) 

BPRr3 

1.0 + (Ag3 byPASs/^S3 CORE^ 



X 


^S3 BYPASS CX)RE 
^R3 BYPASS ^R3 CORE 


- 1.0 
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For core stream rotor 3 


1 + BPRr3 

As3 bypass ^S3 CORE 

1.0 + ^Ag3 bYPASS^^3 core) 

^R3 BYPASS ^R3 CORE 


(A13) 

For core stream stator 3 




^R4 '^S3 CORE 
As3CORE 


(A14) 


Two values of mach number are calculated in each row. The axial mach number is the mach 
number in a direction parallel to the centerline of the compression system. To determine the 
total pressure in each row, the mach number in the direction of the flow must be deter- 
mined. The axial mach number is determined from the isentropic flow relationship 


nO.5 


M 


1.0 + |W„ „/Pso „ A„ MOD21 2 12.0 (i,„ - 1.0)R „/g T„1 - 10 


yc n 


Tn-'O 


(A15) 


The flow direction Mach number (M 2 ) is a function of the axial Mach number and flow angle 
(a). A schematic of the flow angles is shown in Figure A-6. 

The flow angle at the exit of a stator is assumed to be the stator exit metal angle and is a 
constant value. The flow an^e at the exit of a rotor is influenced by the rotational speed 
of the rotor and therefore mu.st be calculated. The method of computation is 


V 


yen 




y (l.O/M^yp jj) + (7n - 1.0)/2.0 8in^ ac n 


(A16) 


and 


iina 


cn 


1.0 

^/1.0 + (U„A^ycn“»0/taiiU 

7n^ 


(A17) 
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Vgg u is the axial flow velocity and ^ is the exit metal angle of the rotor. The relationship 
of these parameters is shown in Figure A-7. 

Equations (A16) and (A17) are solved by substituting equation (A17) into equation (A16) 
which provides a quadratic equation to solve for and subsequent calculation of 

®c n* 

The flow mach number is calculated from the axial mach number and the flow angle as 


^zc n ^yc n^*^ “ 


c n 


(A18) 


The ccnnpressor flow function (0^) contains a total pressure term which must be deter- 
mined along the flow path. The total pressure is a function of the flow mach number and 
b related to the static pressure using the isentropic flow relationship 

yn 


^Tc n “ ^Sc n 


1.0 + 


(7n-10)M",cn 

2.0 


7 „- 1.0 


(A19) 


The flow conditions across the gap are assumed to be constant with one exception. The 
static pressure and the flow angle, and therefore flow mach number and total pressure, 
are modified at the exit of core rotor 3, shown in Figure A-3. The modification is to account 
for the splitter flow separation at low bypass ratio. 


The modifications are implemented through a change in the parameter ^ which b the core 
rotor 3 metal exit angle. The modified angle, Ug b defined as a function of bypass ratio. 
The value of „ is 


^g „ = 75.5277 - 14.3077 BPR (A20) 

with a minimum limit of the actual core rotor 3 metal exit angle (57.5“). 

The axial velocity at the gap exit is recalculated by 


'^yg n "■ ^yc n ^g n^^” ^c n 


The static pressure at the exit of the gap is determined by 


**Sg n “ **Sc n 


2.0 W 


ac n 


‘n g 


V -V 
yc n yg n 


1.0 + A„+| MOD2/A„ 


(A21) 


(A22) 
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MOD2 is the core stator 3 area modification presented previously. 
The flow angle at the gap exit is recalculated as 


VygnWn-Vn/«"Oc«> 

v'1.0 + CVyg„/(U„-Vy^„/tmfl,„)l2 


The change in the axial velocity and the flow angle results in a change in the flow mach 
number which is determined by 


nil'll 

10 - (-r„ - 1.0) „/2.0 y„ RgTTj „ « g „ 


The total pressure at the exit of the gap is calculated from the isentropic flow relationship 


/ (T„-10)m2 \ 

. P 3 ,„ (..0 


Tn-10 


(A25) 
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Burner 


Description 

The TF30-P-3 combustor is composed of eight cans. Air discharged from the high compres- 
ses is divided into a primary and secondary zone. The primary air supports combustion and 
the secondary air provides cooling for the cans and dilution of the fuel-air mixture. 

Simulation 

The burner model used in the engine simulation was separated to simulate the dynamics 
of the hot combustion section and the cool bypass flow. Separation of the burner primary 
and secondary airflow was therefore necessary to simulate this operation. Two parallel air- 
streams were used to represent the primary and secondary zones, with the combustion cal- 
culations performed in the primary airstream. At the discharge plane of the burner, the 
primary and secondary streams were mixed to estaUish the inlet conditions to the turbine. 
The model of the primary zone of the burner accounted for combustion efficiency, heat 
addition, pressure losses, and mixing of the two zones, 'fhe secondary zone is represented 
as a duct, which is described under the section entitled “Ducts”. The burner is shown in 
Figure A-1 where station 4.0 is the burner entrance and station 5.0 is the burner exit. 

Efficiency 

During the combustion process, all of the theoretical temperature rise is not obtained from 
the fuel-air mixture. The actual temperature rise is taken into account through use of an 
efficiency factor which is empirically derived from both rig and engine tests of the burner. 

The efficiency characteristics used in the simulation were applied to the heat addition 
computations of the primary zone. The efficiency was described in curve form as illustrated 
schematically in Figure A-8 where 


t7b = f(f/a, Pt4.0> ^T4.0’ EFF^ 


(A26) 


The burner efficiency constant Kb efF “ 1.7 and has been determined 

experimentally from engine data. 

Heat Addition 

The primary zone heat addition process was simulated in terms of temperature rise. The 
temperature rise in the burner is determined using the enthalpy of formation of a compound. 
The technique used for chemical reactions is to determine the enthalpy of the reacting 
components, fuel and air in the burner, at some reference state (25**C(77*F) and 1 atmos- 
phere pressure). The enthalpy of combustion is the enthalpy change when combustion 
occurs and the state of the elements is changed from the reactants to the final combustion 
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products. Since initially the reactants are not at the reference state, a computation must be 
performed to determine the enthalpy change required to attain the reference state. After 
the combustion has occurred, the enthalpy change effects a temperature change in the final 
combustion products. 

The basic method of implementing the temperature rise is 


^^REACTANTS ^REACTANTS "" ^PRODUCTS 


where Hj^y^CTANTS “ energy required to cool or heat the reactants to the reference 
state. 


“reactants = “ 4.0 PRI • “reference state reactants (A28) 

The energy release due to the combustion of the fuel is determined from the enthalpy of 
combustion for the fuel (T^mBUSTION^* flowrate and the burner efficiency 


QrELEASE " ’lB^f*'COMBUSTION 


(A29) 


The energy release plus the energy required to cool or heat the reactants is used to in- 
crease the temperature of the products of combustion. The enthalpy of the products can 
be determined by 


“ 5.0 PRI = ^RELEASE + “REACTANTS + “REFERENCE STATE PRODUCTS 

The temperature rise is determined from the enthalpy and the pressure specific heat 
coefficient for the products. The pressure specific heat coefficient is a function of tempera- 
ture therefore an integration process must be used to determine the final temperature of 
the products of combustion as shown 


•*5.0 PRI “ **5.0 PRl/^a4.0 PRI ^f^ 


where 


Cp ■ f(TT5 Q pRi) 


r ^T5.0 PRI 

•I CpdT 

^ 25"C (77®) 


(A31) 
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Pressure Drop 

KT.:rrJsr^ S“ ~ 

*T4.(/^STA ND^ 
.o/**STAND / 


APb = Pt4.qKb 


^Bi Wiiivil 15 . 

AaPRI^ 
\ *"t4.o/**I 


(A33) 


The (XMistant Kg is determined experimentally. 

Mixina 

The primary and secondary zones of the burner are mixed at the exit of the burner The 

of «>' n.on»„t»m «kI en«*y. i„ eqw 


mass 


'^aPRI+W3SEc = WaS.o 


(A34) 


momentum 

*a PRI ’PRI + Ps PRI ApRi + W, SEC vjEc + Pj jec *SEC - 
^aS.O '^S.O **S5.0 ^5.0 


(A35) 


energy 


*a PRI 1«B PRI + Wa SEC "b SEC * *»$.o "bS-O 


(A36) 


A further constraint which must be satisfied is that the 
two streams to be mixed must be equal. That is 


sUtic pressures at the exit of the 


**SPRI*PSSEC 


(A37) 
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Turbine 


Description 

The TF30-P-3 engine uses a single stage hi^ turbine to power the high compressor and a 
three stage turbine to power the fan-low compression system. The turbines are shown in 
Figure A-1 with the hi^ turbine between stations 5.0 and 6.0 and the low turbine between 
stati(Mis 6.0 and 6.9. 

Simulation 


The performance of the low and high turbines is determined in the same manner. Individual 
stage characteristics for the turbines were not available, therefore it was necessary to re- 
present each turbine from overall performance data. The turbine model accounted for the 
effects of efficiency, turbine cooling air, and power extraction. 

The low and high turbine computations are the same therefore the following description 
will describe the high turbine performance but is applicable to the low turbine as well. 

Efficiency 

The required turbine shaft power is determined from the compressor requirement, power 
extracted for external uses, and power required for acceleration or deceleration of the rotor. 
The rotor acceleration or deceleration wiU be discussed under the section entitled “Rotor 
Dynamics”. The required power can be summed up as 


®Tu “ ****COMPRESSOR ™*EXTERNAL * ™*ROTOR (^38) 


The specific work requirement from the turbine is a function of the power and the turbine 
gas flow Wju as follows 


Ahxu = 


Ktu H^*Tu 
Wtu 


(A39) 


where KJ^ is the conversion factor for converting power to an energy term. 

The turbine gas flow is determined from the mixed burner exit flow and accounts for the 
turbine cooling air. 

The efficiency computation uses an efficiency map for determining a reference efficiency 
and applies modifications to account for Reynolds effects and an efficiency scale factor. 
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A schematic of the efficiency map is shown in Figure A-9. The map is entered with an 
actual woric term and a corrected speed term. The constants Cl and C2j|, in Figure A-9 
are defined as 


(Ah/N^)REFERENCE 

(Ah/N^) design 


and 


(A40) 


^T5.0^ REFERENCE , ^ ^ ^ 

C2ju = ■ ^ (A41 ) 

(N/^ Tj5 q) design 

The efficiency map is a generalized map which can be used for a variety of engines. The 
reference values in equations (A40) and (A41 ) are a constant value for the generalized map 
and are determined when the map is generated. The design values are associated with the 
particular engine for which the generalized map is being used and allow the generalized map 
to be used for the TF30-P-3 simulation. In effect, the constants scale the efficiency map to 
simulate various turbine sizes. 


The output of the reference efficiency map is an ideal work term which can be represented 
as shown in Figure A-9. The reference efficiency is 


(CljuAhTu/N^) actual 
(CljuAhTu/N^) ideal 


(A42) 


The efficiency correction for Reynolds effects and scaling are applied to the reference 
efficiency to determine the turbine efficiency 


’^Tu \u REF ^’^u REYNOLDS Tu SCALE FACTOR (A43) 
Turbine Pressure Ratio 


The turbine pressure ratio is detennined from the enthalpy entering the turbine and the 
ideal power extracted from the turbine. The enthalpy at the high turbine entrance 
hjuS.o determined from the exit conditions of the burner. The ideal enthalpy at the 
exit of the turbine is determined from the enthalpy at the entrance, the actual enthalpy 
drop in the turbine, and the turbine efficiency. 


(•^u6.(PlDEAL “ •^u5.0*^*^u%u 


(A44) 
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The turbine pressure ratio is determined from the relative pressure at constant entropy at 
stations S.O and 6.0. The relative pressures are determined by 


**RELATIVE 5.0 “ *’^*»Tu5.0' 

(A45) 

**RELATIVE 6.0* ^[i**Tu6.0>IDEAL' 

(A46) 


The turbine pressure ratio is equal to the ratio of the relative pressures or 


PRTu 

Turbine Area 


^TS.O 

**T6.0 


^RELATIVE S.O 
**RELATIVE 6.0 


(A47) 


The turbine area is calculated by 


^Tu 


^^a5.0 VT5.0 \/PPcRITICAl1 


T5.0 *^CRITICAI 


FP 


(A48) 


The critical flow parameter, P^CRITICAL’ “ d*tennined from starulard gu dynamic rela- 
tionships for a choked nozzle. The ratio of critical flow parameter to turbine flow parameter 
is obtained from a flow parameter ratio map. A schematic of the map is presented in 
Figure A-10. The input to the map is a pressure ratio function and a speed function. The 
constants C3 jy and C4J^ are defined as 


C3tu = 


(N/jlyg q) referen ce 
(N/ ^5.0) design 


(A49) 


and 


C4t„ - <‘o- l o/raTu) referen ce 
( 1.0- *-0/PRxu^ DESIGN 

The definitions of the reference and design values for the flow parameter ratio map are the 
same as previously stated for the efficiency map. 
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Exhaust Nozzle 


Description 

. • -ui- f>r»a Hovice which is used to control the ai«flow 
The TF30-P-3 exhaust noale is a vanable "IL „„ „pe„,e either 

through the engine. The noale is shown m Figure A-1 . The noirle op 

choked or unchoked. 

Simulation 

TWO curves are used to deHne the ope-tion of 

curve, shown schematically ^ ’ ** . effective area to the geometric area. A 

discharge coefficient is defined as the ra i A* 12 is also a function of the 

flow parameter ratio curve, shown schematically in Fi^reA 12 

nozzle pressure ratio. The geometric area is determined by 


^ 10.0 


_ W3to.ojT7^ gPcRlTlC^^ 


1.0 


Cd P 


TIO.O 


FP 


CRITICAL, 


(A51) 




PPCRITICAL “ 10.0^ 


(A52) 


Control 

Description 

Simutation 

TTia simulation of the control schedule uses a curve of the fom, shown in Figure A-13. 

The FLA input determina the Wfffg output is multiplied by 

y® thw required ^operate the engine at steady 

Po to obtain Wf output. It the luei no .nrine accelerates. As the engine 

sSte. excess powa b •™'>®Vm'’‘d™ Sli. t^rdrS^sS^^^^ the steady sUte 
accelerates and N2 increases. ** , . .-aine occurs if the Wf/P» is too low 

operating ,LT “o®” 
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limit the rate of acceleration of the engine and a deceleration schedule is used to limit the 
rate of deceleration of the engine. These limits are shown on Figure A- 13. 

Ducts 

Description 

Five sections of the flowpath in the engine are simulated. These are: the fan duct which 
extends from the fan bypass exit to the exhaust nozzle, the core duct which extends from 
the core flitter to the exhaust nozzle, the turbine exhaust duct which extends from the 
low turbine exit to the core splitter, the diffuser duct which extends from the high compres- 
sor exit to the burner entrance, and the secondary zone of the burner (burner bypass). No 
augmented operation was investigated during this program, therefore, only the pressure loss 
characteristics required simulation. 

Simulation 


The pressure loss was calculated as a function of the duct mach number as shown. 


^DUCT 

*"duct 


KjM2 


Th« constant Kq is a value which was determined from experimental data. 


(A53) 
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DYNAMIC EQUATIONS 


This section describes the equations used to simulate dynamic operation of the engine. 

The equations describe the dynamics of rotor speed transients, exhaust stream mixing, un- 
steady compressor aerodynamics, and fluid flow. 

Rotor Dynamics 

Simulation of engine transient operation requites accounting of the effects of rotor dynam- 
ics. Acceleration or deceleration of the rotors exists when the turbine hcHsepower either 
exceeds or is less than the required steady state level. 


The basic formula used to determine the rotor speed changes is 
NI (dN/dt) 


***‘rotor " 




(A54) 


or rearranging 



(A55) 


Dynamfc Mixin g 

Mixed flow turbofan engine simulations can use a mixing routine which forces an axial 
static pressure balance between airstreams at the exhaust mixing plane. Experience at Pratt 
& Whitney Aircraft indicates that this approach is sufficiently accurate for simulation of 
steady state operation and engine power transients. 

Simulation of engine ( ' 'h frequency operation requires a dynamic mixing routine. The 
dynamic mixing routine used in the modified engine simulation assumes that the inertial 
and viscous characteristics of the gas flow induce a lag in the static pressure balance at the 
mixing plane. At the mixing plane, the dynamic mixing routine separates the duct and core 
airstreams. A static pressure imbalance is allowed to occur between the two airstreams dur- 
ing a transient. The pressure imbalance changes the area split at the exit of the two streams 
in a direction which tends to equalize the pressures. The exit area of each aintream and 
flie static pressure imbalance are related by the following equations: 

^bypass * ^bypass initial value GAIN /<fs bypass “ ^Score^*** (A56) 

Acore " Atotal ~ (A57) 

The magnitude of the gain level in the above equation wm emperically set to give frequency 
reqranse characteristics similar to engine test results. 
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Unsteady Compressor Aerodynamics 

The ctnnpressor unsteady aerodynamics account deviations in rotor and stator lift char- 
acteristics which occur during operation in a dynamic flow field. The unsteady efTects are 
computed by applying a first order lag to the quasi steady entropy change across the rotor 
or stator. The quasi steady entropy change is computed fixmi the equation.' 

As = — R£nTRj - RfinPRj (ASS) 


where TRj and PRj arc the total temperature and total pressure ratios defined by the 
rotor or stator characteristics. The quasi steady entropy change is modified by a time 
constant which is a function of the rotor or stator chord length and the axial velocity 
at the stator or rotor exit. 


diord length 
Axial velocity 


(A59) 


The dynamic entropy change is obtained by applying the time constant as shown below. 

As^l— :^Ls' (A60) 


ri + 1 


Equation (ASS) is rewritten to solve for the dynamic pressure ratio 



£n TRy ■ 



(A61) 


where the temperature ratio is held constant at the quasi steady value. 
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Fluid Flow Dynmic Equations 


The general relationships which describe one-dimensional fluid flow phenomena recognize 
three basic properties. These are the mass of the fluid, its momentum and its energy. These 
jMX>perties are conservative in that a change to any of them within the boundaries of a 
system must be the result of quantities passing through the boundaries. In the most basic 
terms, then, the relationships can be expressed as: 


AMass(^V “ in" out 

Momentum Qy = Momentum jn - Momentum out 

AEnergy cv " Energy - Energy out 


(A62) 

(A63) 

(A64) 


where: 


CV refers to a bounded system known as a control volume or finite element; 
in, refers to quantities entering the finite element; out, to quantities leaving the 
finite elr .nent. 


A refers to an incremental change 


The introduction of the independent variable time (t), allows the formulation of equations 
that describe the properties within the finite element at a particular value of time. 


mt 


Masst = Masst_^t 


AMassin At 
At 


At 


(A65) 


= Momentum^ = Momentumt_^t 


A Momentum^ At 
At 


AMomentumout 

At 

AEnergyjjj At AEnergy^^l At 

E* ' Energy* « Energy*^* 

* At At 


(A66) 

(A67) 


As At approaches zero, the ri^t hand sides of these equations are the integrals with 
respect to time of die time derivatives of the incoming and outgoing properties. These 
are expressed as follows: 



(A68) 


P 


(A69) 


At this point, it is recognized that the solutions obtained must be in terms of the parameters 
with which fluid systems are normally identified. These are pressure, temperature, density 
and velocity. First, formulas are presented fw the time derivatives of the three basic proper- 
ties at the boundaries of the system. This will be followed by formulas for the quantities 
of the properties within the finite element at a particular time. 

With reference to Figure A- 14, a representation of the geometry of a typical finite element, 
the time derivatives of the properties, either incoming and outgoing, are given by: 


d m 


p Av = W 

d t 

(A71) 

BP B (mv) vBm mBv 

= = + - vW + 2F 

B t B t Bt Bt 

(A72) 

BE 

= hyW 

Bt * 

(A73) 

where: 


hT * /cp dT + 0.5 pv2 = Jcp dTj 

(A74) 


The mass derivative is the density times the area perpendicular to the flow direction times 
the distance travelled by the fluid per time differential. This latter quantity is the velocity 
of the fluid. The time rate of change of momentum consists of two terms, the first one 
describing the change in momentum due to mass transfer and the second one consisting of 
die summation of forces on the fluid in the direction of the flow. The second term represents 
the application of Newton’s second law to the finite element. The energy derivative is 
expressed in terms of the total enthalpy which is the sum of the specific kinetic energy of 
the fluid and the specific potential, or thermal, energy contained within the fluid. The 
^cific enthalpy is multiplied by the flow to obtain the time rate of change of energy at 
the boundary. 

At this point, the time derivatives of mass, momentum and energy have been identified 
in terms of the parameters usually associated with fluid flow. In particular, these quantities 
have been defined at the boundaries of a finite element or control volume assuming that 
the flow field is one-dimensional. Spatial integrals which define the quantities of mass, 
momentum and energy which lie within the boundaries of a finite element are given by: 

"*= (/pdV = y^AdL (A75) 

** -^AvdL (A76) 

^ v2 + Cv TI pAdL (A77) 
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A basic assumption of our implementation of the finite element technique is that the 
par^eters of the element vary in a linear fashion from one side to the other. These formula 
tor the spatial vanations in the general parameter arc: 


P “ Pin 



L 


(A78) 


where L varies from 0 to AL across the element This formula is applied to the variations in 
pressure, temperature, density, area and velocity across the element. 


Substitution of these equations into the spatial inugrals allow their evaluation yielding the 
tollowing functions for the properties within the finite element. 


n* = f(p. A, AL) 

P =f(p,A.v,AL) 


(A79) 

(A80) 


E =f(p,A,v,T,AL) 


(A81) 


One item must still be defined. That is the sum of the forces on the fluid exerted by its 
surroundings in the flow direction (2F), which appears in the equation for the time rate of 

reference to Figure A-14, it can be seen that the geometry 
m the fimte element plays a role in the assessment of the forces exerted on the fluid. The 
boundanes trough which flow passes are assumed to be perpendicular to the flow direction. 
The walls of the duct form the other boundary of the element. The total force consists of 
two pai^; the force exerted by the fluid outside of the element and the force exerted by 
the walls on the fluid inside the element. The first term is simply the pressure times the 
flow area at the boundary. In the absence of any area variation across the element or friction, 
there are no body forces exerted by the walls in the flow direction. 


The sum of forces on the fluid within the finite element is then: 


2F = Pin Ain - Pout Aout (A82) 

Up to this point, the discussion has been limited to one-dimensional flow in ducts without 
considering the changes to the properties imposed by external mechanisms; in particular 
changes due to the forces on the fluid exerted by the turbo-machinery and the rate of 
energy transfer to the fluid in the form of mechanical work. Denoting these parameters as: 

F - force on the fluid due to the mechanical components 
dW 

- rate of mechanical work pcrfonned on the fluid by the mechanical 
components 
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the spatial integrals of Equations A75 through A77, define the three equations wWch govern 
the fluid flow within a one dimensional finite element These equations can be written as. 


*** “ (2 Pin A|n + 2 Pout ^out ^in ^out ^out ^in^ - 

P = [ (3 A|n + Aout) ^in '^in (^in ^ ^out) ^out '^out (^in ^out) 

(Pin vout + Pout Vin) 1 AL/12 =/vff.^ y.^ - + 

Pin ^in “ Pout ^out (A84) 

E = Ivin^ (Ain (Pin/5 + Pout/20) + Aout (Pin/20 + Pout/30)) + 

Vout^ (Ain (Pin/30 + Pout/20) + Aout (Pin/20 + Pout/5)) + 

Vin Vout(Ain (Pin/10 + Pout/15) + Aout (Pin/*5 +Pout/10)) 1 AL/2 + 

Cv ( (3 Ain + Aout) Pin ^in (Ain 3 Aout) Pout "(^out 


(Ain Aout) (Pin '^'out Pout '(^in)l AL/1 2 
f dW 

j in ^in “ out ^out ■*" 


(A85) 


In evaluating the external force and energy derivative inputs to a finite element, a concept 
known as the actuator disc is employed. In this model, the turbomachinery is represented 
by a series of discs of negligible length across which total pressure and temperature change 
as a function of a specific component characteristic. The effect of these changes on the 
momentum and energy flux at the incoming boundary of the finite element is included 
by utilizing the parameters that result from the component characteristics in the calcu- 
lation of the incoming momentum and energy derivatives. This concept is shown schemati- 
cally in Figure A- 15. The “in” station represents the inlet to the component while, at the 
“in'” station, conditions that would exist at the component exit in steady-state are 
assumed. Momentum and energy derivatives computed as a function of conditions at 
“in'” will include the effects of the turbomachinery components. On this basis, the 
expressions for the time derivatives of momentum and energy given in Equations (A84) 
and (ASS) can be modified in the following maimer: 


)^in* Vin' "*■ Pin* Ain ® ^in Vin ^in Ain P 


(A86) 


hrin'V=hTi„Wi„ + ^ 


(A87) 
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The modifled engine simulation described the compressor characteristics on a row by row 
baas. Dynamically, the compressor operation was described using ten finite elements. The 
finite element configuration required that row performance characteristics be combined 
algebraically to establish the properties associated with the finite element actuator disc. 

SOLUTION TECHNIQUE 

Three ordinary differential equations (ODE’s) which are applicable to a single finite element 
within an arbitrary flow field were defined in the previous section. This set can be applied 
to each of the elements of the system and results in a set of equations which, when solved 
simultaneously, will describe tlie flow field at any point in time, given the time history of 
the field up to that point and its current boundary conditions. The simulated response of 
the system is composed of a series of solutions at the discrete values' Of time required by 
the particular problem being studied. It must commence with a steady state solution which 
provides the required initial values of the properties. 

The transient solution technique applies implicit solution methods directly to the ODE’s. 

A simple trapezoidal numerical integration scheme is applied to the time derivatives of the 
properties (Equation A71-A73). This results in a set of non-linear algebraic equations which 
yield solutions at discrete values of time in an iterative manner. 


The concept referred to above is generally known as a relaxation method. Relaxation 
methods rely on solving a set of implicitly defined relationships in an iterative manner. The 
solution method uses a set of independent variables which will permit computation of 
both the right hand and left hand sides of the set of governing equations (Equations A83 - 
ASS). A solution is obtained when, for a particular set of boundary conditions, the value of 
all independent variables satisfy the equations. 

In applying this method to the set of non-linear algebraic equations which are used to 
model the dynamic behavior of the flow field described by an arbitrary set of finite elements, 
a set of three variables describing conditions at the exit of each element must be chosen. 

Two of these must describe the state of fluid while the third must describe its motion. Pressure 
and temperature, either static or total, are chosen for the first two while mass flow or 
velocity is used for the third. This set of variables is sufficient to explicitly calculate both 
sides of the three equations and, hence, the three error terms associated with them. 

Additional pairs of independent variables and errors are used to implicitly define boundary 
conditions or constraints of the system that cannot be expressed explicitly. These include 
the flow areM in the turbines and at the exhaust nozzle. A sophisticated non-linear equation 
solver, utilizing a modified Newton-Raphson approach, operates on this set of independent 
variables and errors to achieve the required solutions. 


There are several advantages associated with the use of implicit solution techniques. The 
most important of these is the improvement in execution time possible relative to explicit 
numerical integration techniques. In finite element models, the degree of improvement is 
a function of the number of elements used to represent the system. Explicit techniques 
require a time increment which approaches an order of magnitude smaller than the time to 
propagate a signal across the element. Implicit methods, because of their intrinsic computa- 
tional stobiUty, can tolerate any value of time increment. The selection of the proper one 
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is based purely on the required frequency content of the results. Typically, implicitly defined 
models achieve execution times ten times faster than their explicit counterparts for com- 
parable accuracy. 

Another significant, although often overiooked, advantage to implicit solution techniques 
is the ability to easily modify the boundary conditions of the problem. For instance, 
steady-stote solutions are easily obtained by forcing the time derivatives of the fiuid pro- 
perties to zero. This is accomplished by automatically specifying the time derivatives to be 
the error terms. The set of independent variables remains unchanged. In a similar manner, 
the linearized partial derivatives used in the transfer function calculations are obtained by 
holding the fluid properties (Mass, Momentum and Energy) constant at their given base 
values while each of them are perturbated in turn. In this case, the error terms are the 
difference between the given values of the properties and those calculated from the model . 
Non-zero time derivatives result from each perturbation solution and are used to compute 
the required partial derivatives. The importance of this capabUity lies in the fact that the 
same computer program can be used for all of the different requirements without resorting 
to expensive and time-consuming recoding. 
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TF30-P-3 


TURBINE EXHAUST DUCT 



INLET GUIDE VANE BURNER PRIMARY ZONE 


Figure A-l Engine Station and Component Identification 
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Figure A-2 Flow Parameter Curve 








Figure A-4 Rotor and Stator Static Pressure Characteristics 
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Figure A-5 Rotor Total Temperature Characteristics 
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F^re A-6 Stator Exit Velocity Triangle 


Figure A-7 Rotor Exit Velocity Triangle 




F^re AS Burner Efficiency Map 


Figure A-9 Turbine Efficiency Map 
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Figure A-tO Turbine Flow Parameter Map 


Figure A-1 1 Exhaust Nozzle Discharge Coefficient 
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Ftgure A~12 Ejduaut Nozzle Flow Parameter 
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Figure A-13 Engine Control Schedules 



Figure A-14 Finite Element Representation 
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Figure A-I5 Finite Element Representation Showing the Actuator Disc Concept 
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APPENDIX B 

SIMULATION FREQUENCY RESPONSE 

Frequency response characteristics were estabUshed for the engine and compi^on system 
simSlatiOTs fo?the purpose of evaluating simulation dynamic characteristics, methods 

phase and gain characteristics of the simulations, '^e^ were: 

1) linearized transfer functions, 2) step response and 3) discrete frequency 
general, results obtained by analysis of simulation transiente (step 
Ly differ from those obtained from transfer functions. 

confined to the linear range of the simulation, the difference is attnbutable to the nature 
of transfer function and transient solutions. Transfer function results 
expUcit solution of the system of linear equations derived from the simulation The transfer 
?unc«^lL is ex.^ in the sense that it contains no errors associated with numm^ 
integration. Transient results do include the errors associated with numencal integration. 
?^eTSnces in results can be viewed as a measure of how well the transients reflect the 
dynamic characteristics of the simulation transfer functions. 

In the piocess of obtaining the matrices of linear coefficients which were used ^ 
transfer functions, some observations were made which bear discussion. It was 
the transfer function frequency response results were sensitive to program iteraUon toler- 
»c« and the perturbatior. size of the input and state variables. The sensitrvity was Aown 
by large offsets in the frequency response curves, especially at the low “ 

soectrum The offsets were essentially eliminated by perturbing the variables in both a 
Srandliegative direction, instead of just one direction. A possible 
hnproved results is that the errors associated with the iteration tolerances tended to^cel 
with the two directionally opposed perturbations. In addition, the directionally oppose 
perturbation effectively doubled the perturbation size at a fixed level of iteration 

tolerances. 

The linear matrix coefficients derived from two cases using different perturbation sizes 
were examined and revealed that, although the frequency response characteristics were m 
agreement, the values of the coefficients were different. This indicated the relative con- 
sistency of the coefficients was more important than the levek. As a ^ 

frequ^cy effects noted above appeared to be a function of the relative values of the entire 

set of coefficients. 


A study of the effects of iteration tolerance and pertubation size on frequency 
accuracy and program execution time was made using the compression system simulaticm. 

The result, of Uie study ere shown on Figure B-1. The figure *»*VSVh^d eti e “ 
Size and perturbation size in terms of the percentage error in gam at 0. 1 he^ and case ex 
ecution time in units of minutes and seconds. The ratio of perturbatiOT to tolerance size 
U also shown for each case. Three trends in results are shown by the figure. First, me pro- 
gram execution time increases as the iteration tolerances are reduced. 

!^uirement to make more iterative passes to satisfy solution tolerances. Second, the 
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frequency response accuracy improves as the perturbation to tolerance ratio increases. As 
ttie ratio increases, the influence of the errors introduced by iteration tolerances is dimished* 
ed. Third the perturbation size reaches as optimum with respect to frequency response 
accuracy. The perturbation size must be large enou^ to tes(riye the characteristics of the 
system, but small enough to remain within the linear range of the system. 

A final observation is made on the basis of experience gained from computing transfer 
functions for die compression system and engine simulations. The difficulty of ccmiputing 
transfer functions appeared to correlate with the non uniformity of dynamic element sizing. 
Results of the reduced order compression system investigation showed two cases where 
not all of the iterations were completed within the ranges of input tolerances. The dynamic 
configuration for both cases contained elements of non uniform sizing. The computation of 
transfer functions for the modified engine simulation was also found to be difficuit and 
required linearization of die performance characteristics and scaling of tolerances to obtain 
marginally satisfactory results. The engine simulation was again nonniniform in element 
sizing. The non uniformity was a consequence of introducing the modified compression 
system model into die engine simulation. 
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Figure B-1 Summary of Transfer Function Pertubation Size and Tolerance Study 









LIST OF SYMBOLS 


A 

A 

b 

B 

BPR 

c 

C 

C1,C2, C3.C4 

Cd 

Cp 

CpR 

Cv 

^WAC 

db 

DT 

D 

e 

E 

f/a 

F 

F' 

FP 

g 

h 

H 

HP 

Kb 

Kd 

I 

ID 


Constant 
Area, cm^;in^ 

System matrix 
Constant 
Input matrix 

Bypass Ratio, bypass airflow/core airflow 

Constant 

Output matrix 

Constants 

Exhaust nozzle discharge coeflident 
Spedflc heat at constant pressure, J/(kgXK); Btu/(lbm) (*R) 
Compressor map pressure ratio at intersection of in lines 
Spedfic heat at constant volume, J/(kgXK); Btu/(lbm) (®R) 
Compressor map corrected airflow at intersection of m lines 
Decibels, 20 log (output/input) 

Computer simulation time increment, seconds 
Output/input direct coupling matrix 
Base of natural system of logarithnu 
Energy, (kg) (cm^)/sec^; (Ibm) (in^)/sec^ 

Fuel-air ratio 
Forces, N;lbf 
Mechanical forces, N ; Ibf 

Flow parameters, (Kg) (K)*/^/(scc^) (N); (lbm)(*R)*^2/(secKlb0 
Gravitational constant, l(KgXm)/(NXscc2);32.17(lbmXft)/(lbO(8ec^) 
Spedflc enthalpy, J/Kg; Btu/lbm 
Total enthalpy; J/sec; Btu/iec 
Power, W; Horsepower 

Burner pressure loss coeffldent, sec^/Kg^; sec^/Ibm^ 

Duct pressure loss coeflident 

Moment of inertia, (NXmXsec^); (IbfXFtXsec^) 

Inside diameter, cm; in 


ISO 



LIST OF SYMBOLS (Cont’d) 


» 


IGV 

Kr 

KTu 

L 

AL 

m 

m 

M 

n 

N 

N1 

N2 

OD 

ODE 

OPl 

OP2 

P 

P 

P 

PLA 

PR 

PSTAND 

Q 

R 

R3 

s 

5 

S3 

t 


Inlet guide vane 

Compressor modifier coefficient 

Conversion constant, l(NXm)/(secXW); 550(ftKlbO/(sec) (horsepower) 

Horsepower to energy conversion factor, 1 J/sec; 0.707 Btu/sec 

Distance, cm; in 

Element length, cm; in 

Mass, Kg; Ibm 

Compressor map “m” line 

Mach number 

Stage number designation 

Rotational speed, RPM 

Low rotor speed, RPM 

High rotor speed, RPM 

Outside diameter, cm; in 

Ordinary differential equation 

Operating point at N 1 approximately equal to 8600 RPM 

Operating point at Nl/v/^^j^approximately equal to 7600 RPM 
General designation of finite element parameter 
Pressure, N/cm^; Ibf/in^ 

Momentum, (KgKcm)/sec; (ft)(lbm)/sec 
Power lever angle, degrees 
Pressure Ratio 

Standard pressure, 10.1325 N/cm^; 14.696 Ibf/in^ 

Energy, J/sec; Btu/sec 

Gas constant, 287 (NKm)/(KgXK); 53.3 (lbfXft)/(lbmX“R) 

Rotor three 
Entropy, J/g; Btu/lbm 
LaPlace Operator 
Stator three 
Time, seconds 
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LIST OF SYMBOLS (CONT*D) 


T 

TR 

’^STAND 

u 

U 

V 

V 
W 
W 

WAC 

WAE 

WAT 

X 


Temperature, K; 

Temperature ratio 

Standard temperature, 288.2K; SIS.T^R 
Input vector 

Mean rotor speed, cm/sec; in/sec 
Velocity, cm/sec; in/sec 
Volume, cm^; in^ 

Wei^t flow, Kg/sec; Ibm/sec 
Work,J;Btu 

Corrected airflow, Kg/sec; Ibm/sec 
Engine core airflow. Kg/sec; Ibm/sec 
Total airflow, Kg/sec; Ibm/sec 
State vector 


X 


dx/dt 


X 

y 

a 

P 

A 

5 

n 

T 

X 

t 

0 

p 

T 

e 


Independant variable 
Output vector 
Flow angle, degrees 
Rotor metal exit angle, degrees 
Incremental change 

Ratio of total pressure to standard pressure 

Efficiency 

Specific heat ratio 

Temperature function 

3.141S9 

Flow function Kg/(secXcm^); lbm/(aecXin^) 
Pressure function 
Density, Kg/cm^; Ibm/ft^ 

Time constant, sec 

Ratio of total temperature to standard temperature 
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LIST OF SYMBOLS (CONT’D) 


Subscripts 


a 

Airflow 

B 

Burner 

c 

Compressor stator or rotor exit 

CV 

Control volume 

D 

Duct 

DBS 

Design 

EFF 

Efflciency 

f 

Fuel 

g 

Compressor gap exit 

in 

Entrance of finite dynamic element 

in^ 

Exit of actuator disk 

n 

Stage number designation 

out 

Exit of finite dynamic element 

PRI 

Primary 

AP 

Pressure loss 

R 

RotOT 

REF 

Reference 

S 

Static 

STAND 

Standard 

SEC 

Secondary 

T 

Total 

Tu 

Turbine 

t 

Time 

y 

Axial direction 

z 

Flow direction 

1 

Time designation, seconds 

2 

Time designation, seconds 

3 

Time designation, seconds 

2.0 

Engine inlet 

2.1 

Fan first stator passage 

2.3 

Fan core exit 

2.3F 

Fan bypass exit 


UST OF SYMBOLS (CONT’D) 


3.0 

Low compressor exit 

3.12 

Twelfth stage sUtor passage 

4.0 

High compressor exit 

5.0 

Burner exit 

6.0 

High turbine exit 

6.9 

Lo w turbine exit 

7.0 

Turbine diffuser duct exit 

7.0F 

Fan bypass duct exit 

10.0 

Exhaust nozzle 
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